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Abstract.  Errors  often  occur  in  transferring  electronic  data,  ranging  from  sensitive  gov¬ 
ernment  information  to  everyday  bar  codes.  Encoding  information  with  an  error-correcting 
code  can  alleviate  the  problem  of  corrupt  or  lost  data.  In  order  to  not  overburden  comput¬ 
ing  systems,  an  efficient  code  must  be  used  that  will  quickly  encode  and  decode  data  while 
detecting  and  correcting  a  large  number  of  errors.  The  goal  of  this  project  is  to  construct 
and  develop  efficient  codes  using  recent  advances  in  algebraic  geometry,  combinatorics,  and 
commutative  algebra. 

The  mathematics  of  subspace  arrangements  and  simplicial  complexes  lend  themselves 
well  for  applications  to  coding  theory.  A  subspace  arrangement  is  a  finite  collection  of 
subspaces  in  a  vector  space.  A  simplicial  complex  is  an  abstract  generalization  of  a  polygon 
or  Euclidean  solid.  Fortunately,  both  simplicial  complexes  and  subspaces  arrangements  can 
be  described  algebraically  by  a  collection  of  polynomials,  which  can  be  used  to  construct 
a  code.  Then  the  combinatorial  and  geometric  properties  of  subspace  arrangements  and 
simplicial  complexes  can  be  used  to  enumerate  these  efficient  codes. 

Scripts  and  algorithms  were  written  in  the  computer  algebra  systems  Sage  and  Macaulay2 
to  compute  properties  of  the  codes.  The  data  led  to  the  main  results  of  the  project:  formu¬ 
las  for  the  the  length,  dimension,  and  minimum  distance  of  polygon  and  skeletal  simplicial 
complex  evaluation  codes.  Scripts  were  written  that  aided  in  the  construction  of  proofs 
for  these  formulas.  The  formulas  give  favorable  lengths  (short  to  minimize  computation), 
dimensions  (large  to  allow  for  more  codewords),  and  minimum  distances  (large  to  allow 
more  errors  to  be  corrected  and  identified)  of  these  polygon  and  skeletal  simplicial  complex 
evaluation  codes. 

The  last  part  of  the  project  involved  extensions  to  a  cryptosystem  based  on  these  codes. 
A  cryptosystem  deals  with  enciphering  a  message,  which  is  an  algorithmic  process  designed 
to  make  a  sent  message  unreadable  to  an  interloper,  but,  after  another  algorithmic  deci¬ 
phering  process,  is  readable  to  the  intended  receiver  (who,  unlike  the  interloper,  knows  the 
deciphering  key).  Work  was  done  on  extensions  to  the  code-based  McEliece  Cryptosystem, 
which  has  been  demonstrated  to  withstand  theoretical  quantum  computing  attacks  that 
would  render  common  modern  ciphers  useless. 


Keywords 
Coding  Theory 
Cryptology 

Subspace  Arrangements 
Simplicial  Complexes 
McEliece  Cryptosystem 
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1.  Introduction 

As  error-correcting  codes  have  developed  over  the  past  century,  so  have  the  mathematical 
disciplines  of  algebra,  geometry,  and  combinatorics,  such  as  the  development  of  the  properties 
of  a  characteristic  polynomial  describing  a  poset  representing  a  subspace  arrangement  by 
Athanasiadis  in  [1],  The  goal  of  this  project  is  to  construct  and  develop  efficient  codes  using 
recent  advances  in  algebraic  geometry,  combinatorics,  and  commutative  algebra.  In  particular, 
the  mathematics  of  subspace  arrangements  and  simplicial  complexes  lend  themselves  well  for 
applications  to  coding  theory.  A  subspace  arrangement,  A,  is  a  finite  collection  of  subspaces 
in  a  vector  space,  V.  A  simplicial  complex,  A,  is  an  abstract  generalization  of  a  polygon  or 
Euclidean  solid.  Fortunately,  both  simplicial  complexes  and  subspaces  arrangements  can  be 
described  algebraically  by  a  collection  of  polynomials  (namely,  the  ideal  Ij i),  which  can  be 
used  to  construct  a  code.  Let  C(A,  j)  be  the  evaluation  code  for  the  arrangement  A  with 
polynomials  up  to  degree  j. 

A  major  problem  in  linear  coding  theory  is  finding  codes  that  have  a  small  number  of 
digits  (length)  with  a  high  number  codewords  (dimension),  as  well  as  good  error-correction 
properties  (minimum  distance).  In  addition,  for  certain  classes  of  codes,  finding  the  length, 
dimension,  and  minimum  distance  becomes  very  difficult.  Then  aim  of  this  project  is  find 
the  length,  dimension,  and  minimum  distance  of  codes  of  the  form  C(A,j)-  Fortunately,  the 
combinatorial  and  geometric  properties  of  subspace  arrangements  and  simplicial  complexes 
can  be  used  to  find  the  length,  minimum  distance,  and  dimension  of  these  codes. 

The  formula  for  the  length  of  C(A,  j)  for  any  subspace  arrangement  A  can  be  extracted  from 
the  characteristic  polynomial  of  the  arrangement’s  intersection  lattice.  However,  in  general, 
finding  the  dimension  and  minimum  distance  seems  to  be  a  much  more  difficult  task.  In  Section 
5,  we  find  the  dimension  of  C(A,  j)  for  certain  cases  where  A  is  the  arrangement  of  coordinate 
hyperplanes.  We  then  focus  on  the  binary  simplicial  complex  codes  C(Aa,j)-  The  dimension 
for  C(AaA )  is  given  by  the  face  vector  of  A.  Nonetheless,  the  minimum  distance  is  still  very 
difficult  to  determine.  In  Section  6.2,  we  give  an  explicit  formula  for  the  minimum  distance  of 
C (Aai  j)  f°r  all  values  of  j  when  A  is  a  polygon. 

The  most  significant  class  of  codes  we  study  is  the  skeletal  codes  C(Aa,J)  where  A  is 
a  skeleton  (a  simplicial  complex  with  all  possible  faces  of  a  certain  dimension  and  all  lower 
dimensions) .  These  codes  turn  out  to  be  punctured  Reed-Muller  codes,  as  well  as  being  related 
to  Hamming  codes.  The  main  result  is  Theorem  6.22,  which  gives  a  formula  for  the  minimum 
distance  of  C(AaA),  a  skeletal  code  with  j  =  1.  Though  the  minimum  distance  for  j  >  1 
appears  to  be  very  difficult,  we  prove  a  few  cases  of  a  presented  conjecture. 

One  component  of  the  project  involved  the  production  of  computer  scripts  and  algorithms 
in  Sage  and  Macaulay2.  For  example,  a  Macaulay2  algorithm  allowed  the  construction  of  a 
code  (from  which  the  properties  of  the  code  could  be  calculated),  whereas  another  algorithm 
in  Sage  created  a  code  but  only  returned  the  properties  of  the  code.  Other  scripts  helped 
speed  up  lengthy  calculations,  making  proofs  easier  to  see  and  construct. 

The  remainder  of  the  project  was  focused  on  the  McEliece  cryptosystem.  A  cryptosystem 
deals  with  enciphering  (rather  than  encoding)  a  message,  which  is  an  algorithmic  process 
designed  to  make  a  sent  message  unreadable  to  an  interloper,  but,  after  another  algorithmic 
deciphering  process,  is  readable  to  the  intended  receipt  (who,  unlike  the  interloper,  knows 
the  key  to  decipher  the  message).  The  McEliece  cryptosystem  is  built  on  error-correcting 
codes  and  is  suitable  to  codes  with  nice  decoding  algorithms  (see  [6],  [5],  and  [18]).  Since  our 
subspace  arrangement  codes  are  closely  related  to  Hamming  codes  and  Reed-Muller  codes, 
it  is  possible  that  they  also  have  nice  decoding  properties.  An  important  note  here  is  that 
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the  McEliece  cryptosystem  has  been  demonstrated  to  withstand  some  theoretical  quantum 
computing  attacks  that  would  render  common  modern  ciphers  (such  as  RSA)  useless. 

2.  Basic  Definitions 

In  this  section  we  establish  notation  and  the  basic  objects  of  study. 

2.1.  Coding  Theory.  We  begin  by  systematically  defining  all  objects  of  study  for  this  project 
with  a  focus  on  linear  codes  as  the  main  character  of  study.  To  do  this  we  define  all  relevant 
terms. 

A  commutative  ring  is  a  set  of  elements  closed  under  two  operations  (usually  called  addition 
and  multiplication)  that,  under  addition,  forms  an  Abelian  group  (a  set  of  elements  with  an 
operation  that  has  a  unit  and  that  is  associative  and  commutative,  and  inverses  exist)  and, 
under  multiplication,  forms  a  commutative  monoid  (a  set  of  elements  with  an  operation  that 
has  a  unit  as  well  as  being  associate  and  commutative;  inverses  do  not  necessarily  exist).  The 
distributive  property  also  holds  in  a  commutative  ring.  A  field,  is  a  commutative  ring  where 
all  multiplicative  inverses  exist.  A  finite  field  is  a  field  with  finitely  many  elements.  We  call  Fg 
the  finite  field  with  q  elements  and  note  that  q  must  be  a  power  of  a  prime.  One  example  is 
F2,  which  is  defined  as  the  set  {0, 1}  with  addition  and  multiplication  modulo  2.  This  example 
is  of  prime  importance  because  of  its  applicability  to  computers. 

A  vector  space ,  V,  is  a  set  of  vectors,  {iq, . . .  }  closed  under  addition  and  scalar  multiplica¬ 
tion,  with  the  distributive  property  holding.  The  central  example  we  study  is  the  vector  space 
F” .  A  subset  of  a  vector  space  which  is  closed  under  vector  addition  and  scalar  multiplication 
is  a  subspace  of  the  vector  space.  The  minimum  number  of  vectors  required  to  generate  all 
vectors  of  a  subspace  [or  vector  space]  (through  scalar  multiplication  and  vector  addition)  is 
the  dimension  of  the  subspace  [or  vector  space].  For  example,  the  dimension  of  F”  is  n. 

The  aim  of  coding  theory  is  to  develop  efficient  data  transmission  processes  that  will  detect 
and  correct  errors.  A  code  consists  of  a  collection  of  codewords,  each  of  which  are  assigned 
to  a  symbol  (or  group  of  symbols).  A  simple  example  is  the  ISBN  10-digit  number  used  for 
cataloging  books.  The  first  nine  digits  are  used  for  identification,  and  the  last  digit  is  a  linear 
combination  of  the  identification  digits  that  checks  for  errors  when  the  bar  code  is  scanned 
(here,  the  codewords  are  the  10-digit  numbers  assigned  to  each  book). 

In  order  to  introduce  coding  theory  terminology,  consider  the  following  example.  Let  M  = 
{a,b,c,d}  be  a  set  of  information.  A  code  for  M  could  consist  of  the  subspace  of  codewords 
C  =  {(0, 0, 0,  0),  (0, 1, 0, 1),  (1, 0, 1,  0),  (1, 1, 1, 1)},  with  each  codeword  representing  one  element 
in  M.  Note  that  the  codewords  in  this  example  all  have  digits  in  F2.  Formally,  a  code  C  is  a 
subspace  of  a  vector  space  V  =  F”  over  the  finite  field  F9;  in  this  example,  C  is  a  subspace  of 

V  =  F|.  The  length ,  n,  of  a  code  refers  to  how  many  digits  are  in  a  codeword,  which  is  also 

the  dimension  of  the  ambient  vector  space  V ;  in  the  example,  n  =  4.  The  dimension ,  k.  of 
a  code  is  the  dimension  of  the  subspace  C;  in  the  example,  k  =  2.  A  basis  of  vectors  for  the 
subspace  C  provides  for  the  construction  of  the  generating  matrix  of  C ,  which  is  the  matrix 
with  the  basis  vectors  forming  the  rows.  In  the  example,  the  generating  matrix  for  C  is 

T  1  0  1  0 

G  “  [  0  1  0  1 

Finally,  the  minimum  distance ,  d ,  is  the  minimum  number  of  digits  that  must  be  changed  to 
transform  one  codeword  into  another;  in  this  example,  d  =  2.  Minimum  distance  can  also 
be  defined  as  the  minimum  Hamming  weight  among  the  non-zero  codewords  (the  Hamming 
weight  of  a  codeword  is  the  number  of  non-zero  entries  in  it).  Minimum  distance  dictates  the 
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number  of  errors  a  code  can  detect.  In  particular,  a  code  can  detect  up  to  d  —  1  errors  and 
correct  up  to  |_fj  errors  (note  that  [z\  is  the  largest  integer  less  than  or  equal  to  z).  A  code 
with  length  n,  dimension  k ,  and  minimum  distance  d  is  referred  to  as  an  [n,  /c,  d]q  code,  where 
the  digits  are  elements  of  Fg.  The  example  is  a  [4,  2,  2] 2  code. 

2.2.  Hamming  Codes.  An  important  class  of  codes  which  will  relate  to  the  codes  developed 
in  this  project  is  the  class  of  Hamming  codes.  All  binary  Hamming  codes  can  be  constructed 
in  a  fairly  straightforward  process,  which  is  described  in  [9].  For  the  Hamming  code  denoted 
l~Lr,  first  form  a  matrix  with  the  numbers  1  through  2r~1  in  binary  as  the  columns  (so  there 
are  2r  —  1  columns).  Permute  the  matrix  to  form  a  matrix  of  the  form  [A|/r],  where  Ir  is  the 
r  x  r  identity  matrix  (this  matrix,  termed  Hr,  is  called  the  parity-check  matrix  of  the  code). 
Then,  T~Lr  can  be  described  by  the  matrix  [/2>-_i_r|  AT],  where  AT  is  the  transpose  of  A  (the 
columns  of  AT  are  the  rows  of  A ,  with  the  ordering  of  the  rows  becoming  the  ordering  of  the 
columns).  Each  binary  Hamming  code  'Hr  has  the  property  of  being  a  [2r  —  1,  2r  —  1  —  r,  3] 2 
code. 

Example  2.1.  The  construction  of  %4  is  provided.  First,  the  matrix  of  binary  numbers  is 

'000000011111111' 
000111100001111 
011001100110011  ' 
101010101010101 

The  columns  are  then  permuted  to  form 

'110100011111000' 

_  110111110000100 
111001100110010  ■ 
101110101010001 

From  this  matrix,  the  generating  matrix  for  H4  is  formed: 

'100000000001111' 
010000000001110 
001000000000011 
000100000001101 
000010000000101 
n4  =  000001000000110  . 

000000100000111 
000000010001100 
000000001001001 
000000000101010 
000000000011011 

By  inspection,  it  can  be  confirmed  that  H4  is  a  [15, 11,  3]2  code. 

2.3.  Reed-Muller  Codes.  Another  class  of  closely-related  codes  is  the  class  of  Reed-Muller 
codes,  which  are  described  in  detail  in  [14].  A  binary  Reed-Muller  code,  7 Z(j,m),  is  formed 
by  examining  all  polynomials  in  m  variables  up  to  degree  j.  Each  polynomial  corresponds  to 
a  codeword  by  evaluating  that  polynomial  on  all  possible  points  in  F™.  For  each  7 Z(j,m),  a 

[2m,  J2  (™)>  2m~r]i  code  is  produced. 
i= 0 
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Example  2.2.  The  construction  of  the  72.(1,  3)  Reed-Muller  code  will  be  demonstrated.  For 
consistency,  the  following  ordering  of  the  points  in  F3  will  be  used: 


Xi 

£2  £3 

0 

0  0 

0 

0  1 

0 

1  0 

0 

1  1  . 

1 

0  0 

1 

0  1 

1 

1  0 

1 

1  1 

Each  codeword  is  obtained  by  creating  a  binary  string  that  represents  the  evaluation  of  each 
polynomial  on  the  certain  states  of  x\,  x2,  and  £3,  with  the  order  given  by  the  above  table. 
Thus,  the  following  table  can  be  made  to  display  all  the  codewords: 

Polynomial 

Codeword 

0 

(0,0, 0,0, 0,0, 0,0) 

Xi 

(0,0, 0,0, 1,1, 1,1) 

x2 

(0,0, 1,1, 0,0, 1,1) 

X3 

(0,1, 0,1, 0,1, 0,1) 

X\  +  £2 

(0,0, 1,1, 1,1, 0,0) 

X\  +  £3 

(0,1, 0,1, 1,0, 1,0) 

£2  +  £3 

(0,1, 1,0, 0,1, 1,0) 

£1  +  £2  +  £3 

(0,1, 1,0, 1,0, 0,1)  . 

1 

(1,1, 1,1, 1,1, 1,1) 

1  +  £1 

(1,1, 1,1, 0,0, 0,0) 

1  +  £2 

(1,1, 0,0, 1,1, 0,0) 

1  +  £3 

(1,0, 1,0, 1,0, 1,0) 

1  +  £1  +  £2 

(1,1, 0,0, 0,0, 1,1) 

1  +  £1  +  £3 

(1,0, 1,0, 0,1, 0,1) 

1  +  £2  +  £3 

(1,0, 0,1, 1,0, 0,1) 

1  +  £1  +  £2  +  £3 

(1,0, 0,1, 0,1, 1,0) 

Note  that,  since  most  of  the  polynomials  are  linear  combinations  of  one  another,  72.(1, 3)  can 
be  represented  by  a  matrix  of  a  basis  for  the  code  (the  polynomial  is  in  the  first  column): 

1  11111111 
xi  00001111 
£2  00110011' 

£3  0  1  0  1  0  1  0  1 

1 

Note  that  the  code  is  a  [23,  (3),23_1]2  =  [8, 4, 4]2  code. 

2.4.  Subspace  Arrangements.  Let  V  be  a  vector  space  of  dimension  £  over  a  finite  field 
F  of  q  elements.  A  subspace  arrangement  A  =  {Xi, . . . ,  Xt}  in  V  is  a  finite  collection  of 
linear  subspaces  Xi  C  V.  Some  mathematical  objects  related  to  subspace  arrangements  are 
now  introduced  (definitions  for  some  of  them  are  contained  in  the  following  sections).  These 
objects  are  further  discussed  in  [2].  Let  S  =  Fg \x-\_, . . . ,  X(\  be  the  symmetric  algebra  of  the  dual 
vector  space  V*  (the  dual  vector  space  is  the  set  of  all  linear  maps  </>  :  V  — >  Fg).  The  points 
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Figure  1 .  The  coordinate  planes 


of  V  in  A  are  P(A)  =  (J  Xi  =  {p\, . . .  ,pn}-  Additionally,  let  1(A)  =  {/  £  S\f(P(A))  =  0} 
2—1 

be  the  defining  ideal  of  A.  With  these  properties  of  a  subspace  arrangement,  it  is  possible  to 
develop  an  error-correcting  code  from  a  subspace  arrangement. 

3.  Subspace  Arrangement  Codes 

3.1.  Definition  of  C(A,  j).  In  order  to  define  the  codes  we  want  to  study  we  need  to  under¬ 
stand  polynomials  and  their  roots  over  a  finite  field.  A  polynomial  ring,  Fg[®i, . . .  ,xt\,  is  the 
set  of  all  polynomials  over  the  variables  x\ ,xt  with  usual  addition  and  multiplication  and 
with  coefficients  in  Fq.  A  subring  of  a  ring  is  a  subset  of  a  ring  that  is  closed  under  the  ring 
operations.  An  ideal,  I,  is  a  subring  that  is  closed  under  multiplication  by  any  element  in  the 
ring. 

Recall  the  notation  of  subspace  arrangements.  We  can  now  construct  a  code  from  A  by 
evaluating  all  polynomials  of  certain  degrees.  Define  the  evaluation  map  evA  '■  S<j  — >  F™  by 

evA(f)  =  (f(Pl),---,f(Pn)) 

where  S<j  is  the  vector  space  of  all  polynomials  of  less  degree  than  or  equal  to  j  in  S  = 
Fg[xi, . . . ,  xe].  Now  we  can  define  our  main  object  of  study. 

Definition  3.1.  The  image  C(A,j )  =  im(evA)  is  a  linear  subspace  in  F”  that  we  call  a 
subspace  arrangement  code. 

Now  we  examine  a  few  examples. 

Example  3.2.  Let  q  =  2,  l  =  3,  j  =  1,  and  A  be  the  xy-,  xz-,  and  yz-planes,  as  seen  in 
Figure  1.  Then,  1(A)  =  <  X\X2X3  >.  Thus,  V  =  F|  and  S  =  F q[x\X2Xz).  It  follows  that 
P(A)  =  {(0,0,0),  (1,0,0),  (1,1,0),  (1,0,1),  (0,1,0),  (0,1,1),  (0,0,1)},  so  \P(A)\  =  7.  Now,  to 
find  \m(evA  '■  — >  Fj),  we  write  a  basis  for  the  subspace  spanned  by  the  image  of  the 

evaluation  map  as  a  matrix  (the  top  row  designates  a  point  in  P(A),  and  the  first  column 
delineates  the  polynomials  in  S<  1  at  which  each  point  is  evaluated): 


(0,0,0) 

(1,0,0) 

(1,1,0) 

(1,0,1) 

(0,1,0) 

(0,1,1) 

(0,0,1) 

1 

1 

1 

1 

1 

1 

1 

1 

Xi 

0 

1 

1 

1 

0 

0 

0 

X2 

0 

0 

1 

0 

1 

1 

0 

x3 

0 

0 

0 

1 

0 

1 

1 
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Figure  2.  The  xz-  and  yz-planes 

Close  observation  yields  that  the  length  of  the  code  is  7,  the  dimension  is  4,  and  the  minimum 
distance  is  3.  Therefore,  this  code  is  a  [7, 4,  3] 2  code. 

Example  3.3.  Let  q  =  7,  i  =  2,  j  =  2,  and  A  be  the  xz-  and  j/2-planes  as  in  Figure  2.  Then, 
1(A)  =  <  X1X2  >■  Thus,  V  =  Ff  and  S  =  Fyfaqa^]-  It  follows  that  P(A)  =  {(0,0),  (1,0), 
(2,0),  (3,0),  (4,0),  (5,0),  (6,0),  (0,1),  (0,2),  (0,3),  (0,4),  (0,5),  (0,6)},  implying  \P(A)\  =  13. 
Now,  to  find  im(eu.4  :  S< 2  —>  F}3),  determine  its  corresponding  matrix: 


(0,0) 

(1-0) 

(2,0) 

(3,0) 

(4,0) 

(5,0) 

(6,0) 

(0,1) 

(0,2) 

(0,3) 

(0,4) 

(0,5) 

(0,6) 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Xl 

0 

1 

2 

3 

4 

5 

6 

0 

0 

0 

0 

0 

0 

X2 

0 

0 

0 

0 

0 

0 

0 

1 

2 

3 

4 

5 

6 

X? 

0 

1 

4 

2 

2 

4 

1 

0 

0 

0 

0 

0 

0 

xi 

0 

0 

0 

0 

0 

0 

0 

1 

4 

2 

2 

4 

1 

Close  observation  yields  that  the  length  of  the  code  is  13,  the  dimension  is  5,  and  the  minimum 
distance  is  5.  Therefore,  this  code  is  a  [13,  5,  5)7  code. 

Example  3.4.  Note  that  if  A  is  the  entirety  of  V  =  F^  (so  1(A)  =  <  0  >)  and  q  =  2, 
then  P(A)  =  V .  Thus,  the  subspace  arrangement  code  is  equivalent  to  a  Reed-Muller  code. 
Specifically,  the  code  is  C(¥eq,j)  =  7 Z(j,£).  Since  skeletal  codes  (as  discussed  in  Section  6.3) 
are  a  punctured  C(F*,j)  code,  they  are  also  a  punctured  Reed-Muller  code. 

3.2.  Code  Length  and  Characteristic  Polynomial.  Let  A  be  the  subspace  arrangement 
{Xi, . . . ,  A'fc},  where  for  all  i,  X,  is  a  subspace  of  the  coordinate  hyperplanes.  Let  I  a  be  the 
corresponding  ideal  in  F^.  Let  L(A)  consist  of  all  intersections  of  the  subspaces  of  A  (note 
that  the  empty  intersection  is  defined  as  the  entire  vector  space,  V).  Then  L(A)  is  a  lattice 
and  a  poset  by  reverse  inclusion  (a  poset  is  a  partially-ordered  set,  meaning  that  an  order 
relation,  such  as  reverse  inclusion,  exists  on  some  of  the  elements  of  the  set;  a  lattice  is  a  poset 
where  any  two  elements  have  a  unique  least  upper  bound  and  greatest  lower  bound).  Next, 
the  Mobius  function,  fj,,  on  L(A)  is  fi  :  L(A)  — >  Z  defined  recursively  by 

f  H(V)  =  1 
T(X)  =  -  E  n(Y)  ■ 

l  Y<X 

From  this  function,  the  characteristic  polynomial  for  A,  x(*4.,  t),  is  defined  as 

x(A,t)=  MX)tdim(x)). 

XeL(A) 
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Figure  3.  Poset  of  2-Boolean 


In  [1] ,  Christos  Athanasiadis  proved  a  theorem  for  determining  the  number  of  the  set  of  points 
of  A  (in  Fg): 

\P(A)\=q*-X(Aq). 

Since  |P(.4.)|  corresponds  to  the  length  of  a  code,  we  have  the  following  interpretation  of 
Athanasiadis’s  result. 

Corollary  3.5.  If  n  is  the  length  of  a  code  associated  with  the  subspace  arrangement  A,  with 
characteristic  polynomial  x(A,<3'),  then 

n=  |P(A)|  =qe  -x(Aq)- 

The  following  example  illustrates  this  process. 

Example  3.6.  Let  V  =  F§  and  A  be  the  a>axis  and  y-axis.  Then,  1(A)  =  <  xy  >,  with 
a  corresponding  poset  {a,  b,  c,  d},  where  a  <  b,a  <  c,a  <  d,b  <  d ,  and  c  <  d.  For  an 
easier  understanding  of  this  process,  the  poset  can  be  interpreted  to  have  a  correspondence 
to  R2  (and  the  2-Boolean,  which  is  the  Boolean  arrangement  [discussed  in  Section  5]  in  R2) 
as  follows:  a  =  V  =  R2,  b  =y-axis,  c  =x-axis,  and  d  =origin.  Figure  3  is  a  picture  of  the 
Hasse  diagram  of  this  poset.  It  follows  that  L(A)  =  {a,  6,  c,  d}  =  {V,  b,  c,  d}.  Using  the 
Mobius  function,  n(V)  =  1,  p(b)  =  —ply)  =  — 1  (a  <  b ),  p(c)  =  —p(V)  =  — 1  (a  <  c), 
and  p(d)  =  —(p(c)  +  p(b)  +  p(V))  =  —(—1  —  1  +  1)  =  1  (c  <  d,  b  <  d,  a  <  d).  Thus,  the 
characteristic  polynomial  is 

X(A,  t)  =  [1  *  t2]  +  [-1  *  t 1]  +  [-1  *  t1]  +  [1  *  f°]  =  t2  -  2t  +  1 

(each  bracketed  addend  corresponds  to  the  evaluation  of  a,  6,  c,  and  d,  respectively).  Finally, 
evaluating  using  Athanasiadis’s  theorem  gives 

n  =  \A\  =qe-  X(A,  q)  =  32  -  x(A  3)  =  32  -  [32  -  2  *  3  +  1]  =  5. 

Working  out  P(A)  by  hand  confirms  the  answer: 

P(A)  =  {(0,0),  (1,0),  (2,0),  (0,1),  (0,2)}  =►  \P(A)\  =  \A\  =  5. 
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3 


Figure  4.  3-simplex 
4.  Coordinate  Arrangement  Codes 

We  want  to  study  the  dimension  and  minimum  distance  of  C(A,j)  for  any  subspace  ar¬ 
rangement  A,  but  doing  so  is  difficult.  For  the  rest  of  the  paper,  focus  will  be  placed  on 
subspace  arrangements  that  are  coordinate  arrangements,  which  are  described  in  this  section. 

4.1.  Simplicial  Complexes.  In  order  to  study  subspace  arrangements  that  are  intersections 
of  coordinate  hyperplanes,  we  first  define  and  build  a  theory  for  what  are  called  simplicial 
complexes.  We  follow  the  standard  formulation  of  the  correspondence  between  coordinate 
arrangements  and  simplicial  complexes  written  in  [2]  by  Bjorner.  For  convenience,  let  [fc]  = 
{1, . . . ,  k}  be  the  set  of  numbers  1  through  k.  A  simplicial  complex ,  A,  is  a  set  of  subsets  of 
[k]  such  that 

(1)  if  cr  £  A  and  t  is  a  subset  of  a,  then  r  £  A,  and 
(2)  if  a;  £  [k],  then  {a-}  £  A 

where  [k]  are  the  vertices  of  A  and  cr  is  a  face  of  A.  A  fc-simplex  is  a  simplicial  complex  that 
contains  all  subsets  of  [fc].  Geometrically,  we  can  view  the  1-simplex  as  a  point,  the  2-simplex 
as  a  line  segment,  the  3-simplex  as  a  triangle,  the  4-simplex  as  a  tetrahedron,  and  such  forth. 

Example  4.1.  Let  k  —  3  and  A  =  {{1},  {2},  {3},  {1,  2},  {1,  3},  {2, 3},  {1, 2, 3}}.  We  can  think 
of  A  geometrically  as  a  filled-in  triangle  by  associating  a  face  of  size  3  to  the  3-simplex.  In 
fact,  this  simplicial  complex  is  exactly  the  3-simplex  as  pictured  in  Figure  4. 

Example  4.2.  Let  k  =  3  and  A  =  {{1},  {2},  {3},  {1,  2},  {2, 3},  {1,  3}}.  Here  the  simplicial 
complex  is  geometrically  represented  by  a  triangle  with  the  inside  missing,  which  is  pictured 
in  Figure  5. 

The  combinatorial  setting  allows  one  to  easily  study  higher  dimensional  geometric  objects. 
In  order  to  use  simplicial  complexes  to  produce  codes,  we  associate  an  ideal  to  a  simplicial 
complex.  For  a  simplicial  complex  A  with  vertices  [t\ ,  the  vertices  correspond  to  the  variables 
of  the  polynomial  ring  Fg[a:i, . . . ,  X(\,  and  the  ideal  is  I  a  =  {xa  \  cr  ^  A}. 

4.2.  Definition  of  C(Aa,  j)-  The  process  of  correlating  a  subspace  arrangement  with  a  sim¬ 
plicial  complex,  as  described  in  [2],  is  as  follows.  For  {bi,...,b„}  (a  basis  of  F”)  and  each 
subset  a  =  {*i,...,is}  C  [n],  let  Ka  =  span{b(1 , . . . ,  bis }.  The  subspace  arrangement 
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Figure  5.  Empty  triangle 


Figure  6.  Simplicial  complex  corresponding  to  the  x-axis  and  yz-plane 

Aa  =  {Ka\(7  £  A}  is  thus  determined  by  the  simplicial  complex  A.  With  this  correlation 
in  mind,  it  is  now  possible  to  define  a  coordinate  arrangement  code. 

Definition  4.3.  Let  A  be  a  simplicial  complex  with  corresponding  subspace  arrangement  „4a- 
Then  the  coordinate  arrangement  code  corresponding  to  A  is  the  subspace  arrangement  code 

C(AaJ). 

The  construction  of  error-correcting  codes  from  simplicial  complexes,  as  well  as  the  correla¬ 
tion  between  subspace  arrangements  and  simplicial  complexes,  can  be  demonstrated  through 
examples. 

Example  4.4.  Let  A  =  {{1},{2,3}}  C  2I3!.  The  subspace  arrangement  then  consists  of 
the  span  of  the  first  variable  joined  with  the  span  of  the  second  and  third  variables  (that 
is,  <  x  >  U  <  y,  z  >).  This  subspace  arrangement  thus  consists  of  the  rr-axis  together  with 
the  j/z-plane.  The  simplicial  complex  and  subspace  arrangement  are  depicted  in  Figure  6. 
From  the  construction  described  above  and  recalling  the  conventions  discussed  in  Section  3.1, 
it  is  clear  that  £  =  3.  Let  q  =  2  and  j  =  1  (note  that  q  and  j  are  not  contingent  on  the 
construction  of  A),  so  V  =  Fij  and  S  =  F2  [*1^20:3].  It  follows  that  I  a  =  <  X\X2,XiX^  >■ 
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Figure  7.  Tetrahedron 


Hence,  P(Aa)  =  {(0,0,0),  (1,0,0),  (0,1,0),  (0,0,1),  (0,1,1)},  implying  |P(^A)|  =  5.  The 
matrix  corresponding  to  im(ex  :  S< i  — >  F|)  is  given  below. 


(0,0,0) 

(1,0,0) 

(0,1,0) 

(0,0,1) 

(0,1,1) 

1 

1 

1 

1 

1 

1 

£1 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

1 

1 

X2X3 

0 

0 

0 

0 

1 

Thus,  the  length  is  5,  the  dimension  is  5,  and  the  minimum  distance  is  1,  so  the  code  is  a 
[5,  5, 1]2  code. 


Example  4.5.  Let  A  represent  the  simplicial  complex  interpreted  geometrically  as  a  tetrahe¬ 
dron  with  vertices  1,2,3,  and  4  and  all  1-dimensional  faces  (the  lines  connecting  each  vertex), 
as  well  as  the  2-dimensional  face  (a  plane)  {2,3,4}.  Thus,  A  lacks  three  “planes”  and  is  not 
“filled  in”  (since  {1,2, 3, 4}  ^  A).  The  simplicial  complex  can  be  interpreted  geometrically  as 
in  Figure  7. 

From  the  construction  described  above  and  recalling  the  conventions  discussed  in  Section 
3.1,  it  is  clear  that  i  =  4.  Let  q  =  2  and  j  =  1,  so  V  =  F}  and  S  =  F2  [£12:2X3X4].  It  follows 
that  J_4a  =  <  X1X2X3, X1X2X4, X1X3X4  >.  Hence,  P(Aa)  =  {(0,0, 0,0),  (0,0, 0,1),  (0,0, 1,0), 
(0,1, 0,0),  (1,0, 0,0),  (0,0, 1,1),  (0,1, 0,1),  (0,1, 1,0),  (1,0, 0,1),  (1,0, 1,0),  (0,1, 0,0),  (0,1, 1,1)},  im¬ 
plying  |.P(.4a)|  =  12.  The  matrix  corresponding  to  im(ex  :  S<i  — >  Fl,2)  is  given  below. 


(0000) 

(0001) 

(0010) 

(0100) 

(1000) 

(0011) 

(0101) 

(0110) 

(1001) 

(1010) 

(0100) 

(0111) 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Xl 

0 

0 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

X2 

0 

0 

0 

1 

0 

0 

1 

1 

0 

0 

1 

1 

X3 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

0 

1 

X4 

0 

1 

0 

0 

0 

1 

1 

0 

1 

0 

0 

1 

Thus,  the  length  is  12,  the  dimension  is  5,  and  the  minimum  distance  is  3,  so  the  code  is  a 
[12, 5, 3] 2  code. 


4.3.  Upper  Bound  on  Dimension.  In  order  to  find  an  upper  bound  on  the  dimension  of 
a  code,  the  Hilbert  function,  H(M,i )  will  be  introduced.  The  Hilbert  function  deals  with  a 
graded  S'-module,  M.  H{M,i)  is  defined  as  the  dimension  of  Mt ,  the  ith  graded  component 
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of  M,  so  H(M,i )  =dim(Mj).  Furthermore,  if  M  =  C[xi, . . .  ,xe],  then 


H(M,  i) 


i  +  e-i 

£-1 


1 


This  formula  results  from  the  realization  of  (given  how  many  variables  are  in  the  polynomial 
ring)  how  many  different  ways  there  are  of  choosing  a  term  of  a  certain  degree.  The  formula 
follows  from  a  theorem  in  discrete  mathematics  [20],  where  the  number  of  selections  with 
repetition  of  r  objects  chosen  from  s  types  of  objects  is  (r+*_1)  (the  r  objects  are  each  degree 
of  a  term  in  a  polynomial,  and  the  s  objects  are  each  variable).  For  example,  if  there  are  three 
variables,  X\ ,  X2,  and  X3,  and  it  is  sought  how  many  monomials  are  of  degree  2,  then  there 
are  three  slots  to  put  a  total  of  2  degrees,  so  there  are  (3+2_1)  =  (2)  =  ®  different  monomials 
of  degree  2  formed  from  three  variables  ({xf,  x|,  x§,  XiX2,XiX3,  X2X3}).  This  formula  can  be 
applied  in  a  more  general  sense  to  a  simplicial  complex  A. 

The  Hilbert  series  is  obtained  from  the  Hilbert  function.  The  Hilbert  series  of  the  graded 

OO 

module  M  is  defined  as  HS(M,t )  =  H(M,i)t*,  meaning  that  the  Hilbert  series  is  a  power 

»= 0 

series  of  t  such  that  the  terms  of  degree  i  refer  to  the  dimension  of  the  ith  graded  component 
of  M.  For  a  simplicial  complex  A,  vector  space  V  =  and  homogeneous  ideal  /a,  let 
F[A]  =  S/Ia  be  the  Stanley-Reisner  ring.  Since  Ia  is  homogeneous,  F[A]  is  graded.  The  face 
vector  of  a  simplicial  complex  is  defined  by  /*  =  |{cr  €  A  :  |oj  =  i}|  for  0  <  i.  The  face  vector 
refers  to  the  number  of  faces  of  dimension  i  —  1  in  A  (/o  is  defined  as  fo  =  1).  The  Hilbert 
series  of  the  Stanley-Reisner  ring  (as  discussed  in  [17])  is 


00  D-\- 1 


HS(F[A\,t)  =  fo  +  EE 


i= 1  m= 1 


i—  1 

TO  —  1 


J mL  > 


where  D  refers  to  the  dimension  of  A.  As  a  corollary  to  Stanley’s  work  in  [17],  an  upper  bound 
for  dimension  exists. 


Corollary  4.6.  For  a  simplicial  complex  A  with  corresponding  subspace  arrangement  A  a,  an 
upper  bound  for  k  =  dim(C(AAj  j))  exists: 


j  D+l 


k  <  fo  +  EE 


i— 1  m= 1 


i  —  1 
in  —  1 


fm- 


Using  the  combinatorial  formula  discussed  at  the  beginning  of  this  section,  the  upper  bound 
can  also  be  described  as 


An  example  illustrates  this  upper  bound. 


Example  4.7.  For  V  =  F|,  let  A  correspond  to  two  connected  points,  so  A  a  is  a  plane  (and 
the  entirety  of  V),  so  /a  =  <  0  >.  Thus,  fo  =  1,  fi  =  2,  (for  each  of  the  two  0-dinrensional 
points),  and  f^  =  1  (for  the  single  edge).  An  upper  bound  for  the  dimension  of  C{Aa,  1)  can 
be  found.  Using  the  above  formula  for  the  Hilbert  series, 

i  d+i  /  . _ .j  \ 

HS(¥[A],t)  =  fo  +  eel: > 

i—  1  m—  1  '  ' 
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i  1+1 


-  fo + 


i= 1  m—1 


l—  1 

m—l 


fo  +  Yl 


i  —  1 


h 


As  a  check,  note 
to  1:  {1, cci, a:2}- 


=  fo¬ 


il 


h 


—  (1)(1)  +  (1)  (2)  +  (0)(1)  =  3. 

that  there  are  exactly  3  polynomials  in  F2[xiX2]  of  degree  less  than  or  equal 
Thus,  the  upper  bound  is  3. 


5.  Boolean  Arrangement  Codes 

In  this  section,  the  Boolean  arrangement  (consisting  of  all  coordinate  hyperplanes  for  a 
vector  space)  is  used  to  produce  ideals  and  codes.  Notation  concerning  Boolean  arrangement 
codes  follows. 


Definition  5.1.  A  Boolean  arrangement  code,  BC is  a  code  generated  by  the  coor¬ 
dinate  hyperplane  arrangement  and  the  ideal  I  =  <  X\  -  ■  ■  x^>. 

This  section  is  concerned  with  results  dealing  with  the  dimension  of  Boolean  arrangement 
codes  in  certain  circumstances. 

Proposition  5.2.  The  dimension  of  BC(£,1)  is  k  =  t  +  1. 

Proof.  Previous  work  allowed  the  extension  of  a  theorem  by  Richard  Stanley  to  be  applied  to 
this  class  of  error-correcting  codes,  establishing  an  upper  bound  for  dimension: 


The  dimension  of  t  +  1  refers  to  the  number  of  linearly  independent  row  vectors  produced 
by  the  matrix  construction  of  the  evaluation  map  ev^\  for  generating  the  code.  Each  row 
corresponds  to  a  polynomial.  The  polynomials  are  1,  x\, . . . ,  xg,.  Thus,  it  suffices  to  show  that 
the  row  vectors  v0,  Vi, . . . ,  produced  by  the  respective  polynomials  1,  x\, . . . , xg  are  linearly 
independent.  Define  ej  as  the  point  (0,  0, . . . ,  0, 1,  0, . . . ,  0)  in  such  that  all  components  in 
e.j  are  0  except  for  the  ?'th  position.  Furthermore,  define  eg  as  the  origin,  (0,  0, . . . ,  0, . . . ,  0). 
Note  that  regardless  of  the  selection  of  q,  e*  for  i  =  0, . . . ,  i  evaluates  to  0  on  I,  so  e*  €  P(A) 
for  *  =  0, . . . ,  t.  For  simplicity  (but  without  loss  of  generality),  let  eo, . . . ,  correspond  to  the 
first  i  +  1  columns  in  the  generating  matrix  G: 


eo 

ei 

e2 

e^-i 

ec  ■■■ 

1 

1 

1 

1 

1 

1  ••• 

Xi 

0 

1 

0 

0  ••• 

X2 

0 

0 

1 

0 

0  ••• 

. .  1 

0 

0 

0 

1 

0  ••• 

Xi 

0 

0 

0 

1  •  •  • 

G  is  clearly  upper  triangular  and  hence  (by  the  properties  of  an  upper  triangular  matrix)  has 
linearly  independent  rows.  □ 
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Remark  5.3.  For  j  >  1,  the  process  becomes  increasingly  difficult.  Some  results  for  specific 
selections  of  £,  q ,  and  j  follow.  For  results  concerning  the  case  q  =  2,  see  Section  6.3. 

Proposition  5.4.  The  dimension  of  BC(2,2)  whenever  q  >  2  is  k  =  5. 

Proof.  Since  q  >  2,  F2  contains  at  least  3  elements,  0,1,  and  cj.  In  the  matrix  construction, 
let  the  first  5  columns  correspond  to  the  respective  points  (0,  0),  (1,  0),  (0, 1),  (cj,  0),  and  (0,  w), 
and  let  the  rows  correspond  to  the  respective  polynomials  l,Xi,X2,x2,  and  x\.  Thus,  the 
resulting  generating  matrix,  G,  is  as  follows: 


(0,0) 

(1,0) 

(0,1) 

(w,0) 

(0,w)-.. 

1 

1 

1 

1 

1 

1... 

Xi 

0 

1 

0 

CJ 

0... 

X2 

0 

0 

1 

0 

CJ  . . . 

xf 

0 

1 

0 

CJ2 

0... 

xf 

0 

0 

1 

0 

w2... 

We  show  the  rows  are  linearly  independent  by  computing  a  maximal  minor  and  showing  it  is  not 
equal  to  zero.  The  minor  we  examine  is  the  minor  with  columns  1  through  5.  Computing  the 
cofactor  expansion  down  the  first  column  gives  that  the  determinant  is  equal  to  the  determinant 
of  the  following  matrix: 

cj  0 
0  to 
cj2  0 
0  cj2 

The  determinant  of  G'  (and  the  maximal  minor)  is  oj2(oj  —  l)2  ^  0  since  Fg  is  a  field  and  thus 
an  integral  domain,  meaning  it  has  no  zero  divisors  (note  that  w  /  0,  w  /  1).  □ 

Proposition  5.5.  The  dimension  of  BC(3,2)  whenever  q  >  2  is  k  =  10. 

Proof.  First  note  that  10  corresponds  to  the  upper  bound: 


G'  = 


1  0 
0  1 
1  0 
0  1 


to  +  £  — 

TO 


=  10. 


The  upper  bound  corresponds  to  the  polynomials  1,  X\,  X2,  X3,  x2,  x2,  x2,  X1X2,  X1X3,  X2X3. 

Thus,  carefully  selecting  the  points  in  P{A)  for  the  columns  (which  can  be  ordered  in  any 
way),  the  left  part  of  the  generating  matrix  G  is 


1 

xi 


x3 

X1X2 

XiX3 

X2X3 


O  (1,0,0)  (0,1,0)  (0,0,1)  0,0,0)  (0,  u>,  0)  (0,0,  u)  (1,  u>,  0)  (1,0,  cj)  (0,l,u). 
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Then  subtract  rows  1,  2,  and  3  from  rows  4,  5,  and  6,  respectively.  The  resulting  10  x  10 
matrix  to  the  left  is  upper  triangular  with  u2  -  w  on  the  diagonal  of  columns  4,  5,  and  6: 


(0,  0,  0) 

(1,0,0) 

(0,1,0) 

(0,0,1) 

(a,,  0,0) 

(0,  a;,0) 

(0,0,  u) 

(1.W.0) 

(1,0,  u) 

(0,  1,  w)  .  .  . 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1  .  .  . 

XI 

0 

1 

0 

0 

U) 

0 

0 

1 

1 

0  .  .  . 

X2 

0 

0 

1 

0 

0 

U) 

0 

U) 

0 

1 . . . 

x3 

0 

0 

0 

1 

0 

0 

U! 

0 

U) 

U)  .  .  . 

-v.2 

xl 

xh 

x3 

0 

0 

0 

0 

a,2  -  a; 

0 

0 

1 

1 

0  .  .  . 

0 

0 

0 

0 

0 

a;2  —  u> 

0 

a;2 

0 

1 . . . 

0 

0 

0 

0 

0 

0 

2 

U)  —  U) 

0 

a,2 

u;2... 

x\X2 

0 

0 

0 

0 

0 

0 

0 

U) 

0 

0  .  .  . 

x\X% 

0 

0 

0 

0 

0 

0 

0 

0 

OJ 

0  .  .  . 

X2X3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

U)  .  .  . 

Since  q  >  2,  these  entries  (as  well  as  all  other  diagonal  entries)  are  not  zero,  so  the  determinant 
of  the  matrix  (and  the  maximal  minor)  is  non-zero,  directly  implying  that  the  rows  are  linearly 
independent.  □ 


Proposition  5.6.  The  dimension  of  BC(£,2)  whenever  q  >  2  and  £  >  3  is  k  =  +  t  +  1. 


Proof.  First  note  that  k 


P+e 

2 


+  £  +  1  corresponds  to  the  upper  bound: 


m  +  £  — 
m 


(l  +  l)!  (*+l)(l) 

[t  —  1) !  (2)  +  +  2 


1. 


Notice  that  since  this  number  is  the  upper  bound,  the  proposed  dimension  corresponds  to 
all  polynomials  in  Fg[xi, . . .  ,xg\  of  degree  0,  1,  and  2.  Thus,  it  is  sufficient  to  show  that  all 

+  £  +  1  rows  in  the  generating  matrix,  G,  are  linearly  independent.  Since  q  >  2  there 
exists  ueF,  such  that  w  ^  0,1.  We  choose  a  specific  order  on  the  points  of  the  arrangement 
so  that  the  matrix  G  can  be  viewed  as 


O 

o 

o 

o' 

o 

T 

o 

o 

(lwO  ...  0) 

o 

o 

(0  .  .  .  01a;)  .  .  . 

1 

1 

1 

1 

1 

1 

1 

1  .  .  . 

xi 

0 

1 

UJ 

0 

1 

1 

0  .  .  . 

x£ 

0 

0 

0 

U) 

0 

a ; 

U)  .  .  . 

'T2 

X1 

0 

1 

a;2 

0 

1 

1 

0  .  .  . 

X 2 

0 

0 

0 

u2 

0 

a;2  .’  .’  . 

XIX2 

0 

0 

0 

0 

u 

0 

0  .  .  . 

X\X£ 

0 

0 

0 

0 

0 

u 

0  .  .  . 

X£_\  x  £ 

0 

0 

0 

0 

0 

0 

U)  .  .  . 

Subtract  rows  1, ...  ,£  from  rows  £+ 1, . . . ,  2£,  respectively.  The  resulting  leftmost  square  ma¬ 
trix  is  upper  triangular  with  non-zero  entries  on  the  diagonal,  so  it  has  a  non- zero  determinant. 
Thus,  the  maximal  minor  is  non-zero,  implying  that  the  rows  are  linearly  independent.  □ 


Proposition  5.7.  The  dimension  of  BC(4,3)  whenever  q  >  3  is  k  =  35. 


Proof.  Of  first  note  is  the  upper  bound  for  these  parameters: 


3 


On.  4-  £.  —  P 


3 


On.  4-  33 


m— 0 


m= 0 


=  1  +  4  +  10  +  20  =  35. 
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Thus,  in  order  to  show  that  the  dimension  is  35,  it  suffices  to  show  that  each  row  vector  in 
the  matrix  construction  of  the  code  (of  which  there  are  35  for  each  polynomial  of  degree  less 
than  or  equal  to  3)  is  linearly  independent.  Furthermore,  careful  ordering  of  the  columns  of 
the  matrix  construction  will  lend  themselves  to  a  clear  determination  of  linear  independence. 
Let  u  €  F9\{0, 1}  such  that  u  ^  w2  and  oj  ^  — 1.  Thus,  the  matrix  construction  can  be  as 
follows: 


G  = 


Mi  *  * 

0  M2  * 

0  0  m3 


Mi,  M2,  and  M3  are  all  square  sub-matrices  (of  different  dimensions)  along  the  diagonal  of 
G  (thus,  0  and  *  represent  square  sub-matrices  elsewhere  in  the  matrix,  where  0  denotes  all 
entries  in  the  sub-matrix  are  0).  It  will  be  demonstrated  that  Mi,  M2,  and  M3  are  all  able 
to  be  made  to  be  upper  triangular,  meaning  that  G  can  be  made  to  be  upper  triangular  and 
implying  that  all  of  the  rows  are  linearly  independent. 

First,  Mi  contains  the  rows  corresponding  to  the  polynomials  {  1,  Xi,  X2,  x3,  X4,  xf,  x\,  x3, 
x\,  x\,  a;3,  x3,  x\  }  and  the  columns  corresponding  to  the  points  {  (0,0, 0,0),  (1,0, 0,0),  (0,1, 0,0), 

(0,0, 1,0),  (o', 0,0,1)  (w, 0,0,0),  (0,w,0,0)  (0,0, w,0)  (0,0, 0,w),  (w2, 0,0,0),  (0,w2,0,0),  (0,0,w2,0), 

(0,0,0,w2)}.  Observe  that  Mi  is  as  follows: 


Mi 


1111 
0  10  0 
0  0  10 
0  0  0  1 
0  0  0  0 
0  10  0 
0  0  10 
0  0  0  1 
0  0  0  0 
0  10  0 
0  0  10 
0  0  0  1 
0  0  0  0 


1  1  1 
0  to  0 
0  0  w 

0  0  0 
10  0 
0  w2  0 
0  0  w2 

0  0  0 
10  0 
0  w3  0 
0  0  w3 

0  0  0 
10  0 


1  1  1 

0  0  w2 

0  0  0 

u>  0  0 

0  U!  0 

0  0  cc4 

0  0  0 

w2  0  0 

0  w2  0 

0  0  w6 

0  0  0 

w3  0  0 

0  w3  0 


1  1  1 
0  0  0 
w2  0  0 

0  w2  0 
0  0  w2 

0  0  0 
w4  0  0 

0  CO4  0 

0  0  w4 

0  0  0 
w6  0  0 

0  w6  0 
0  0  w6 


In  order  to  make  Mi  upper  triangular,  it  suffices  to  use  row- reduction  techniques  on  M\.  To 
describe  the  row-reduction  techniques  being  used,  some  new  notation  will  be  introduced.  Let 
the  rows  1-4  of  the  above  matrix  be  collective  called  A,  the  next  four  B,  and  the  last  four  C 
(the  0th  row  does  not  need  to  change  during  the  operations  and  thus  is  not  included).  Let 
operations  cA  +  dB  be  defined  as  the  multiplication  of  each  row  in  A  and  B  by  c  and  d, 
respectively,  and  then  the  component-wise  addition  of  A  and  B  (so  the  first  row  in  A  is  added 
to  the  first  row  in  B;  this  process  is  analogous  to  creating  a  new  vector  space).  Thus,  (through 
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row- 

reduction)  replace 

B  by  B 

-  A,  and  replace  C  by  C  —  B  — 

v(B- 

-C): 

"11111 

1 

1 

1 

1 

1 

1 

l 

1 

0  10  0  0 

L J 

0 

0 

0 

0 

0 

0 

0  0  10  0 

0 

OJ 

0 

0 

0 

a;2 

0 

0 

0  0  0  1  0 

0 

0 

OJ 

0 

0 

0 

oj2 

0 

0  0  0  0  1 

0 

0 

0 

OJ 

0 

0 

0 

oj2 

0  0  0  0  0 

2 

UJ  —UJ 

0 

0 

0 

4  2 

OJ  —OJ 

0 

0 

0 

0  0  0  0  0 

0 

OJ2  —OJ 

0 

0 

0 

oj4—oj2 

0 

0 

0  0  0  0  0 

0 

0 

OJ2  —OJ 

0 

0 

0 

4  2 

OJ  —OJ 

0 

0  0  0  0  0 

0 

0 

0 

cu2- 

oj  0 

0 

0 

oj4—oj2 

0  0  0  0  0 

0 

0 

0 

0 

OJ&  —  u;5  —  OJ4  +OJ3 

0 

0 

0 

0  0  0  0  0 

0 

0 

0 

0 

0 

oj6-oj5-cj4+oj3 

0 

0 

0  0  0  0  0 

0 

0 

0 

0 

0 

0  oj6- 

-OJ5  —OJ4 

+ud 

0 

-0  0  0  0  0 

0 

0 

0 

0 

0 

0 

0 

a,6- 

_w*_w4+w3 

Note  that  this  is  an  upper  triangular  matrix,  so  the  determinant  of  it  is  the  product  of  the 
diagonal  entries: 


1 4(w2  —  w)4(w6  —  w5  —  u4  +  CO3)4 


=  W4(  U)  —  l)4^w3((u j  —  l)(w  +  l)(w 


4 


=  w16(w  —  l)12(w  +  l)4. 


Thus,  the  determinant  is  0  if  and  only  if  u>  =  0, 1,  —1,  which  is  not  the  case.  Therefore,  M\ 
can  be  made  to  be  upper  triangular  with  a  non-zero  determinant. 

Now,  examine  M2,  which  contains  the  rows  corresponding  to  the  polynomials  {  X\X2,  X\X 3, 

X\X^,  X2X3,  X2X4 ,  X3X4,  x\x2,  2^2:3,  2:12:4,  X2X3,  x\x4,  x\x4,  X\x\,  X\x\,  X\x\,  2:22:§,  X2X4, 


X3X4  }  and  the  columns  corresponding  to  the  points  {  (1,1, 0,0),  (1,0, 1,0),  (1,0, 0,1),  (0,1, 1,0), 


(0,1, 0,1),  (0,0, 1,1),  (l,w,0,0),  (1,0, w,0),  (1,0, 0,w),  (0,1, w,0), 
(w,0,w2,0),  (w,0,0,w2),  (0,w,w2,0),  (0,w,0,u;2),  (0,0, w,w2)}: 


(0,1, 0,w),  (0,0, l,w),  (w,w2,0,0), 


M2 


1  0  0  0  0  0  w 
0  1  0  0  0  0  0 
0  0  1  0  0  0  0 
0  0  0  1  0  0  0 
0  0  0  0  1  0  0 
0  0  0  0  0  1  0 
1  0  0  0  0  0  OJ 
0  1  0  0  0  0  0 
0  0  1  0  0  0  0 
0  0  0  1  0  0  0 
0  0  0  0  1  0  0 
0  0  0  0  0  1  0 
1  0  0  0  0  0  w2 
0  1  0  0  0  0  0 
0  0  1  0  0  0  0 
0  0  0  1  0  0  0 
0  0  0  0  1  0  0 
0  0  0  0  0  1  0 


0  0  0  0  0 

w  0  0  0  0 

0  w  0  0  0 

0  0  w  0  0 

0  0  0  w  0 

0  0  0  0  OJ 

0  0  0  0  0 

w  0  0  0  0 

0  w  0  0  0 

0  0  w  0  0 

0  0  0  w  0 

0  0  0  0  w 

0  0  0  0  0 

w2  0  0  0  0 

0  w2  0  0  0 

0  0  w2  0  0 

0  0  0  w2  0 

0  0  0  0  w2 


w3  0  0  0 

0  w3  0  0 

0  0  w3  0 

0  0  0  w3 

0  0  0  0 

0  0  0  0 

w4  0  0  0 

0  w4  0  0 

0  0  w4  0 

0  0  0  w4 

0  0  0  0 

0  0  0  0 

w5  0  0  0 

0  w5  0  0 

0  0  w5  0 

0  0  0  w5 

0  0  0  0 

0  0  0  0 


0  0 
0  0 
0  0 
0  0 
w3  0 
0  w3 
0  0 
0  0 
0  0 
0  0 
w4  0 
0  w4 
0  0 
0  0 
0  0 
0  0 
OJ5  0 
0  w5 


In  order  to  make  M2  upper  triangular,  it  suffices  to  use  row-reduction  techniques  on  M2.  Let  D, 
E,  F  represent  (respectively)  the  first  6,  next  6,  and  last  6  rows.  To  make  M2  upper-triangular, 
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switch  E  and  F. 

Then, 

replace  E 

by  E 

- 

D  and  F  by 

F-D: 

100000 

L 0 

0 

0 

0 

0 

0 

a;3 

0 

0 

0 

0 

0 

010000 

0 

to 

0 

0 

0 

0 

0 

a;3 

0 

0 

0 

0 

001000 

0 

0 

to 

0 

0 

0 

0 

0 

CJ3 

0 

0 

0 

000100 

0 

0 

0 

to 

0 

0 

0 

0 

0 

to3 

0 

0 

0  0  0  0  1  0 

0 

0 

0 

0 

to 

0 

0 

0 

0 

0 

a,3 

0 

0  0  0  0  0  1 

0 

0 

0 

0 

0 

to 

0 

0 

0 

0 

0 

CO3 

000000 

2 

CO  —to 

0 

0 

0 

0 

0 

5  3 

to  —to 

0 

0 

0 

0 

0 

000000 

0  to' 

2  —to 

0 

0 

0 

0 

0 

u;5— a,3 

0 

0 

0 

0 

000000 

0 

0 

2 

to  —to 

0 

0 

0 

0 

0  to 

5 —co3 

0 

0 

0 

000000 

0 

0 

0  to 

2  —CO 

0 

0 

0 

0 

0 

to5— to3 

0 

0 

000000 

0 

0 

0 

0  to' 

2 _ 

to  0 

0 

0 

0 

0 

5  3 

to  —co 

0 

000000 

0 

0 

0 

0 

0 

to2— to 

0 

0 

0 

0 

0 

to5 —co3 

000000 

0 

0 

0 

0 

0 

0 

to4— to3 

0 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

0 

0 

0 

M-a,3 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

0 

0 

0 

0  to 

4—co3 

0 

0 

0 

000000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

w4-a;3 

0 

0 

000000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

to4— CO3 

0 

.0  0  0  0  0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a;4-a;3 

The  determinant  is  w24(w  —  l)12,  which  is  nonzero  because  w  ^  0, 1.  Therefore,  Mi  can  be 
made  to  be  upper  triangular  with  a  non-zero  determinant. 

Lastly,  M3  can  be  examined,  which  is  as  follows: 


(1,1, 1,0) 

(1,1, 0,1) 

(1,0, 1,1) 

(0,1, 1,1) 

X1X2X3 

1 

0 

0 

0 

X1X2X4 

0 

1 

0 

0 

X1X3X4 

0 

0 

1 

0 

X2X3X4 

0 

0 

0 

1 

Clearly,  A f3  is  upper  triangular  with  a  non-zero  determinant.  Thus,  G  can  be  expressed  as  an 
upper  triangular  matrix  with  the  determinant  being  non-zero.  Therefore,  all  35  row  vectors 
are  linearly  independent,  so  k  =  35.  □ 

Proposition  5.8.  For  BC(£,j)  with  q  0  (l  >  j),  k  >  £j  +  1. 

Proof.  First,  some  definitions  are  necessary.  Note  that  in  all  of  the  above  matrices  in  the  rows 
which  correspond  to  the  monomials,  for  each  degree  of  each  monomial,  a  diagonal  submatrix 
was  formed.  For  example,  examine  the  following: 


(1,0,0) 

(0,1,0) 

(0,0,1) 

(w,0,0) 

(0,07,0) 

(0,0,  oj) 

Xi 

1 

0 

0 

U! 

0 

0 

X2 

0 

1 

0 

0 

to 

0 

X3 

0 

0 

1 

0 

0 

U) 

xi 

1 

0 

0 

w2 

0 

0 

xi 

0 

1 

0 

0 

cu2 

0 

x\ 

0 

0 

1 

0 

0 

w2 

In  order  to  allow  for  arbitrary  t  to  be  included  in  the  argument,  consider  the  equivalent 
abbreviated  matrix  for  this  matrix  (which  is  not  contingent  upon  the  choice  of  £): 


(i) 

M 

Xi 

1  UJ  ■ 

1 

w2 

Each  entry  in  this  abbreviated  matrix  corresponds  to  an  £  x  £  block  of  the  larger  matrix  of 
the  form  (entry  element)(/^),  where  Ie  is  the  £  x  £  identity  matrix.  If  the  determinant  of  this 
abbreviated  matrix  is  shown  to  be  non-zero  for  certain  q1  then  the  determinant  of  the  larger 
matrix  is  non-zero. 
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Now  we  recall  how  a  determinant  is  calculated.  An  entry  in  the  first  row  of  the  matrix  is 
multiplied  by  1  or  —1  and  then  multiplied  by  the  determinant  of  its  cofactor.  Repeating  this 
process  until  all  cofactor  determinants  are  calculated  reveals  that  this  process  is  analogous  to 
a  permutation.  Thus,  for  a  matrix  A ,  det(A)  =  Y(~ l)aAa,  where  A„  is  a  permutation  and 
a  is  dependent  upon  whether  Aa  is  even  or  odd.  Examine  the  matrix  G  that  consists  of  the 
following  tj  +  1  rows  of  the  matrix  that  generates  the  code: 


(1) 

M 

(ca2) 

(W3) 

A4) 

(W5) 

i 

1 

l 

1 

1 

1 

1 

Xi 

1 

U) 

u2 

w3 

w4 

A5) 

A 

1 

(jj2 

CO4 

w6 

w8 

A10) 

G  = 

A 

1 

w3 

w6 

w9 

w12 

A15) 

A 

1 

w4 

w8 

w12 

w16 

A20) 

A 

1 

w5 

w10 

w15 

w20 

A25) 

a  peri 

nut  at 

ion  a,  for  the  mth  row  (the  initial 

constant  polynomial,  is  termed  the  Oth  row),  a  column  is  uniquely  identified  (which  is  identified 
with  a  term  of  degree  im).  It  follows  then  that  the  total  degree  of  the  term  in  the  maximal 

j 

minor  of  G  for  the  permutation  a  is  Y  (m)(*m)-  The  term  of  highest  degree  (treating  ui  as 

m= 0 

3 

a  variable)  is  sought.  The  main  diagonal  has  the  property  that  ( im )  =  m,  so  Y  (m)(im)  = 

m—0 


Y  m2  =  Al+MM+i)^  Thus,  if  the  main  diagonal  provides  the  highest  degree  term  in  the 

m—0 

j 

maximal  minor,  it  must  be  demonstrated  that  Y  (TO)(*m)  <  where  it  is  assumed 

m—0 

that  s  ^  is  for  at  least  one  such  s. 

Proof  by  induction  will  be  used.  The  base  case  is  j  =  1  (the  j  =  0  case  is  a  lxl  matrix,  so 
the  conditions  of  the  inequality  cannot  be  fulfilled).  The  matrix  for  the  j  =  1  case  is 


(1) 

(w) 

1 

1 

1  . 

Xi 

1  UJ 

The  main  diagonal  term  has  degree  1,  which  is  collaborated  by  the  summation  formula: 

j 

i(i+i)(2+i)  _  6  _  -p  The  other  diagonal  (which  is  the  only  option  for  Y  (m)(*m)  such 

m=0 

3 

that  s  ^  is  for  at  least  one  such  s)  has  degree  Y  (m)(im)  =  (0)(1)  +  (1)(0)  =  0.  Thus,  the 

m—0 

inequality  holds  for  the  base  case. 

j 

Now,  assume  Y  <  j^+1g2j+1)  ■  We  must  show 

m—0 


t+i 

m—0 


(j  +  l)((j  +  1)  +  l)(2(j  +  1)  +  1) 

6 


with  s  ^  is  for  at  least  one  such  s.  Without  the  loss  of  generality,  assume  that  s  is  the  maximal 
index  such  that  is  ^  s.  Thus,  it  =  s  for  some  t  <  s,  and  ir  =  r  for  all  r  >  s. 
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m— 0 
t- 1 


m=t-\- 1 
s — 1 


3  + 1 

Then  ( m)(im )  is  as  follows: 

m=0 

j+l  s  j+1 

E  (m)(im)  =  J2  ( m)(im )  +  0)+)  +  55  (TO)(*m) 

m— 0  m— 0  m=s+l 

m^t 

t—  i  s  j+i 

=  51  (m)(*m)  +  E  (TO)(*m)  +  0)0)  +  E  ^ 

m=s+l 

__  i+i 

=  55(m)(iTO)  +  55(m  +  l)(im+i)  +  Q)p)  +  y  ?n2 

m— 0  m—t  m— s+1 

i-1  s — 1  s —  1  j+ 1 

=  55  (m)(*m)  +  55 (m)(*™+i)  +  55 im+i  +  0)0)  +  55  m2- 

m— 0  m—t  m—t  m=s+l 

Let  {*o,  ■  •  • ,  *t-i,  *t+i,  ■  •  ■ ,  is}  =  {ace  ■  •  • ,  as-i}  C  [s  —  1]  to  combine  the  first  two  summations: 

j  +  1  S— 1  S— 1  j+l 

E++  ^m)  —  E  (■ m)(am )  +  E  ^m+1  H”  (<)(*)+  E  m2. 

m—0  m—0  m=t  m=s+l 

Use  the  induction  hypothesis  on  the  first  term: 

V~V  a/"  \  „  0  ~  00)(2p  -  1)  +  1)  •  i  \  i  E 

2^  (to)  Om)  <  - g -  +  2^  *™+l  +  0)0)  +  2^ 

m—0  m—t  m=s+ 1 

Summation  formulas  can  be  used  to  simplify  the  expression: 
i  ;  i 


to2. 


E+)( 


tm )  < ^ 


m=0 


0-  1)0)(2(g  - 1)  + 1)  0-  i)s  _  0  -  iQ 

6  2  2 


-  ts 


+ 


0  +  l)((i  +  1)  +  l)(2(j  +  1)  +  1)  00  +  l)(2s  +  1) 


6  6 

Thus,  it  suffices  to  show 

P~  l)P)(2p~  1)  +  1)  |  Q~l)s  (t-l)t 


ts 


6  2  2 
(J  +  1)(0  +  1)  +  l)(2(j  +  1)  +  1)  00  +  l)(2s  +  1) 


< 


6  6 

0  +  i)(0  +  i)  +  i)(2p  +  i)  +  i) 


6 


which  is  analogous  to 

p  —  l)p)(2(s  —  1)  +  1)  p  -  l)s  {t—l)t 


6 


+ 


+  ts-  00  +  lp  +  l)  <Q 


and  to 


(sp  +  l)(2s  +  1)  ( (s  —  l)p)(2(s  —  1)  +  1)  |  p  —  l)s  (t  —  l)t 

6 


6  2  2 
Algebraic  manipulation  of  the  left-hand  side  of  the  last  inequality  follows: 

0p  +  l)(2s  +  l)  /  (s  —  l)p)(2(s  —  1)  +  1)  p  —  l)s  (t  —  l)t 

6  V  6  +  2  2 


+  ts  >0. 


ts 
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=  \  Us(s  +  l)(2s  +  1)  —  Us  —  l)s(2s  —  1)  +  3(s  —  l)s  —  3 (t  —  1  )t  +  6 ts))) 

6 

=  -(2s3  +  3s2  +  s  —  (2s3  —  2s  —  3 12  +  3t,  +  6  st))  =  -(3s2  +  3s  +  3 12  —  3 1  —  6  st) 

6  6 

=  ^(s2  +  S  +  t2  -  t  -  2 st)  =  i((s  -  t)(s  -  t  +  1)) 

Since  t  <  s,  s  —  t  >  0,  and  s  —  t  +  1  >  0.  Thus, 

(s(s  +  l)(2s  +  1)  _  / (s-  l)(s)(2(s  -  1)  +  1)  {s-  l)s  _  {t  -  1  )t  \ 

6  \  6  2  2  ) 

=  7>((S-0(S^+1))  >0: 

j  +  1 

directly  implying  Y2  ( m)(*TO )  <  lJ+1)((J+b+1)(2(.?+i)+i) ;  completing  the  proof  that  the  term  of 

771— 0 

maximal  degree  in  the  determinant  of  the  abbreviated  form  of  G  is  obtained  by  multiplication 
along  the  main  diagonal,  which  evaluates  to  .  However,  this  expression  is  not  the 

term  of  highest  degree  in  G,  since  each  term  in  the  abbreviated  form  of  G  describes  an  t  x  t 
submatrix.  Since  this  i  x  t  submatrix  is  the  identity  matrix  multiplied  by  a  constant,  the  only 
way  to  keep  a  term  in  the  determinant  of  G  non-zero  is  by  multiplying  along  the  diagonal  of 
the  submatrix.  Thus,  for  the  submatrix  corresponding  to  a/,  instead  of  contributing  a  term  of 
degree  i  to  a  term  in  the  determinant  using  the  submatrix,  the  submatrix  contributes  a  term  of 
degree  (i)(i)  to  the  term.  Thus,  the  term  of  highest  degree  in  the  determinant  is  gU)U+1)(2J+1)  _ 
It  follows  from  the  existence  of  this  term  that  the  determinant  is  also  not  zero.  By  an  extension 
of  the  Fundamental  Theorem  of  Algebra  to  finite  fields,  there  can  be  as  many  as  £U)U+1)(2j+1) 
zeros  in  the  expression  describing  the  determinant  of  G.  Let  q  >  gU)U+1)(.2J+1)  _  Thus,  there 
are  more  elements  in  the  finite  field  than  zeros  in  the  determinant.  Let  ui  be  an  element  of  the 
finite  field  which  is  not  a  zero  of  the  determinant,  directly  implying  that  the  determinant  of 
G  is  not  zero.  Thus,  all  of  the  row  vectors  of  G  are  linearly  independent.  Therefore,  as  long 
as  q  is  sufficiently  large,  k  >  £(j)  +  1.  □ 

6.  Binary  Simplicial  Complex  Codes 

Error-correcting  codes  are  typically  binary  codes  because  of  the  prominence  of  the  binary 
system  in  computer  operations.  Thus,  in  this  section,  only  binary  codes  ( q  =  2)  are  examined. 
These  codes  derive  from  different  classes  of  simplicial  complexes.  Starting  with  a  proposition 
concerning  the  dimension  for  any  binary  simplicial  complex  code,  the  section  then  moves  onto 
corollaries  concerning  specific  cases  of  these  codes.  An  examination  of  polygons  and  then 
skeletons  (defined  later)  follows.  The  section  ends  in  a  significant  result  regarding  minimum 
distance  for  skeletal  codes. 

6.1.  Dimension.  An  important  part  of  the  examination  of  the  dimension  of  binary  simplicial 
complex  codes  is  the  concept  of  a  face  vector,  /),  in  the  simplicial  complex. 

Definition  6.1.  A  face  vector ,  /),  represents  the  number  of  faces  of  dimension  i—  1  for  a  given 
simplicial  complex.  That  is,  /,;  =  | {cr  e  A  :  \a\  =  i}\  for  0  <  i.  The  trivial  case,  /o,  is  defined 
as  f0  =  1. 

3 

Proposition  6.2.  For  binary  simplicial  complex  codes,  k  =  fi- 

i= 0 
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Proof.  This  proposition  is  analogous  to  the  claim  that  the  evaluation  of  the  Hilbert  function 
upon  a  Stanley-Reisner  ring  based  upon  a  simplicial  complex  (F(A);  the  ring  is  assumed  to  be 
over  Fj  for  some  £)  is  H(SR(A),m)  =  fm  (the  dimension  is  merely  the  sum  for  these  Hilbert 
function  evaluations  from  0  to  j,  which  constitutes  the  upper  bound  previously  determined). 
The  main  reason  that  these  formulas  hold  is  that  in  a  polynomial  ring  over  ¥2,  xn  and  x  (for 
n  a  positive  integer)  evaluate  the  same  over  F2.  Examine  the  formation  of  an  ideal  based 
upon  a  simplicial  complex:  the  ideal  is  based  upon  faces  that  are  not  a  part  of  the  simplicial 
complex.  Thus,  if  a  point  i  and  a  line  segment  jk  are  not  a  part  of  the  simplicial  complex, 
then  x,  and  XjXk  (respectively)  are  contained  in  the  ideal,  meaning  that  the  polynomials  Xi 
and  XjXk  do  not  contribute  unique  row  vectors  in  the  matrix  describing  the  code  (these  row 
vectors  would  be  the  zero  vector  since  the  evaluation  of  the  polynomial  at  each  point  that 
represents  a  column  would  be  zero  since  each  such  point  is  in  P(Aa))-  The  argument  holds 
true  for  cases  for  higher  dimensional  faces.  Thus,  it  readily  becomes  clear  that  the  maximum 
number  of  linearly  independent  row  vectors  in  this  case  directly  corresponds  to  how  many  faces 
of  each  dimension  there  are  in  the  simplicial  complex.  Furthermore,  this  maximum  number 
is  the  exact  number  because  analogously  to  the  previous  proofs  dealing  with  the  Boolean 
arrangement,  points  in  the  ideal  can  be  carefully  chosen  to  create  an  upper  triangular  matrix. 
Therefore,  the  total  number  of  face  vectors  of  dimension  m  that  are  in  the  simplicial  complex 
is  equal  to  the  evaluation  of  the  Hilbert  function  at  m.  Thus,  the  dimension  of  the  code  is  the 
summation  of  the  Hilbert  function  from  m  =  0  to  m  =  j,  which  is  the  summation  of  the  face 

j 

vectors  from  f0  through  fj:  k  =  Y  ,/).  □ 

*= 0 

Corollary  6.3.  Let  M  denote  the  dimension  of  the  minimum  dimensional  non-face  in  a 

j 

simplicial  complex  A.  Then,  if  j  <  M  —  1,  the  dimension  of  C(Aa,  j)  is  k  =  Y  (  ) . 

m— 0 

Proof.  First  note  that  in  F2[aq, . . . ,  X(\,  for  n  a  positive  integer,  the  polynomials  xn  and  x 
evaluate  the  same  on  0  and  1.  Thus,  in  the  matrix  construction  of  the  code,  squarefree 
monomials  need  only  to  be  considered  (any  further  ones  will  result  in  a  row  vector  being 
identical  to  the  row  vector  corresponding  to  the  squarefree  monomial).  Thus,  there  are  (;) 
monomials  with  unique  evaluation  over  F2  of  degree  i.  Since  j  <  M  —  1,  all  of  these  monomials 
with  unique  evaluation  are  in  the  Stanley-Reisner  Ring  (the  faces  corresponding  to  these 
monomials  are  in  the  simplicial  complex).  For  I  =  {*i, . . . ,  it},  let  e/  =  correspond  to 

the  point  in  Fj  with  Is  in  the  i\,. . .  ,it  positions  and  Os  elsewhere.  For  example,  if  /  =  {0,  2, 3}, 
then  e/  =  eo,2,3  =  (1, 0, 1, 1, 0, . . . ).  Let  Xj  be  defined  as  x\  •  ■  ■  Xt  (note  that  the  origin, 
(0,  ...,0)  and  the  constant  polynomial  are  defined  as  e {}  and  Xr\,  respectively).  Thus,  for 

I  7^  {},  Xj(ej)  =  {  J  /  ^  J  '  Now,  construct  G  such  that  the  polynomials  are  normally 

ordered  (from  lowest  to  highest  degree,  in  lexicographic  order  within  the  same  degree).  Order 
the  columns  such  that  if  A/  is  in  the  jth  row,  then  e/  is  in  the  jth  column  (note  that  A  =  {}): 


G  = 


XA 

Xb 

Xc 

XD 


ga 

1 

0 

0 

0 


1 

1 

0 

0 


ec 

1 

* 

1 

0 


eD  ■  ■  ■ 
1 

* 

1 
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Thus,  with  this  careful  ordering,  an  upper  triangular  matrix  is  formed,  thereby  implying  that 
each  of  the  row  vectors  are  linearly  independent.  Since  this  method  of  ordering  covers  all 
monomials  in  the  Stanley- Reisner  Ring,  all  row  vectors  in  G  are  linearly  independent,  so  the 

j 

dimension  of  the  code  is  equal  to  the  number  of  polynomials  up  to  degree  j:  k  =  (m) .  □ 

m— 0 

Corollary  6.4.  If  the  subspace  arrangement  is  the  Boolean  arrangement  (I  =  <  x i  ■  ■  ■  X(  >), 

j 

then  the  length  is  n  =  2e  —  1  and  the  dimension  is  k  =  ^  (  ) . 

m= 0 

Proof.  The  formula  for  the  length  of  the  code  is  trivial:  there  are  a  total  of  2^  points  in  F|, 
and  the  only  point  not  in  the  ideal  is  the  point  with  each  component  being  1,  (1, 1,1,...,  1). 
As  for  dimension,  the  previous  result  can  be  used.  Since  the  Boolean  arrangement  is  used, 

M  =  i.  Clearly,  j  <  i  -  1,  so  k  =  J2  („)•  □ 

m= 0 

Remark  6.5.  While  a  convenient  formula  for  dimension  has  been  fairly  easily  determined,  the 
case  for  minimum  distance  is  much  more  difficult.  In  order  to  obtain  some  results  on  minimum 
distance,  specific  cases  (such  as  letting  j  =  1)  are  examined  in  later  portions  of  this  paper. 

6.2.  Polygons.  A  specific  class  of  simplicial  complexes  are  now  examined;  namely,  the  sim- 
plicial  complex  is  a  polygon. 

Definition  6.6.  Let  Pm  be  the  simplicial  complex  that  is  an  m-gon.  Thus,  Pm  =  {{1},  {2}, 
...,  {?n},  {1,2},  {2,3},  . ..,  {m-  l,m},  {l,m}}. 

Proposition  6.7.  The  characteristic  polynomial  for  a  code  based  upon  Pm  is  x(-Apm,t)  = 
tm  —  mt2  +  mt  —  1 . 

Proof.  Recall  that  L(Apm)  is  a  lattice  that  refers  to  intersections  of  all  subspaces  of  the 
coordinate  hyperplanes  corresponding  to  the  simplicial  complex  A.  The  Hasse  Diagram  for 
L(Apm)  is  depicted  in  Figure  8.  Clearly,  at  the  bottom  of  the  hierarchy  is  Pm  itself,  which 
corresponds  to  the  vector  space  F"'  of  dimension  m.  The  level  of  atoms  consists  of  the  edges 
in  Pm  (which  has  dimension  2),  whose  intersections  from  the  collection  of  vertices  at  the  third 
level  (which  has  dimension  1).  At  the  top  of  the  diagram  is  the  empty  set  (with  dimension 
0),  which  is  the  only  element  in  the  intersection  of  two  vertices.  The  next  step  in  determining 
the  characteristic  polynomial  is  evaluating  the  Mobius  function  upon  the  components  of  the 
diagram  (in  Figure  8,  the  evaluation  of  the  Mobius  function  on  each  component  is  in  red).  By 
definition,  the  Mobius  function  evaluates  to  1  on  Pm.  The  atoms  of  L(Apm)  are  then  each  —1 
(which  is  the  negative  of  the  evaluation  of  the  Mobius  function  on  Pm ).  Each  edge  has  two 
vertices  ordered  below  it,  along  with  Pm,  so  the  Mobius  function  evaluates  to  — (2(— 1)  +  1)  =  1. 
For  the  empty  set,  there  are  m  vertices  (with  value  1)  and  m  edges  (with  value  —1),  along 
with  Pm  (with  value  1),  so  the  Mobius  function  evaluates  to  —  ((m)(l)  +  (m)(— 1)  +  1)  =  — 1. 
Thus,  the  characteristic  polynomial  can  be  formed  by  combing  the  dimension  of  each  level  of 
the  diagram  along  with  evaluation  of  the  Mobius  function:  x(  Apm ,  t)  =  tm  —  mt 2  +  mt  —  1 .  □ 

Corollary  6.8.  The  length  of  C(Apm,j)  is  n  —  2 m+  1. 

Proof.  The  formula  for  length,  when  q  =  2  is 

n  =  2m-  x(MPm,  2)  =  2m  -  ( 2m  -  2 2m  +  2m  -  1) 

=  2m  -  2m  +  4m  -  2m  +  1  =  2m  +  1. 


□ 
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Figure  8.  Hasse  diagram  of  L(Apm ) 


Proposition  6.9.  The  dimension  of  C(Apm,  1)  is  m  +  1,  and  the  minimum  distance  is  3. 


Proof.  Recall  that  k  —  J2  fi ■  Thus,  k  =  J2  fi  =  fo  +  fi  =  1  +  to  =  m  +  1. 


In  F™,  let  e,  =  (0, . . . ,  0, 1,  0, . . . ,  0),  where  the  1  is  in  the  ith  position,  where  i  ranges  from 
0  to  to  —  1  (so  the  initial  position  is  the  0th  position).  For  consistency,  the  origin  is  defined 
as  e{}.  Let  em/j)  then  be  defined  as  e*  +  er  For  Pm ,  Apm  corresponds  to  a  collection  of  to. 
coordinate  planes  in  F ™  which  are  spanned  by  adjacent  coordinate  axis  (that  is,  the  planes 
spanned  by  Xq  and  xi;  Xi  and  X2',  ■  ■  ■ ;  xm-2  and  xm_i;  and  X\  and  xm_i).  The  points  in  the 
subspace  arrangement  are  thus 


p{Apm)  =  1J  Xi  =  {e{}}  U  {a}  U  {e(i)(i+i)} 


X£*Ap7 


for  0  <  i  <  to  —  1.  Note  that  *  +  1  is  addition  modulo  to.  For  simplicity,  let  the  polynomials 
in  the  Stanley-Reisner  ring  up  to  degree  j  =  1  be  labeled  as  1,  Xq,  x\,  . . . ,  xm-2,  and  xm_i. 
For  p  £  Apm, 


1  Tp  6i,  l)(i)}  6(I)(I+1) 

0  else 


Thus,  since  each  non-constant  polynomial  evaluates  to  1  for  only  3  points  in  P(Apm),  the 
minimum  Hamming  weight  among  the  codewords  that  form  the  basis  of  the  code  is  3.  Observe 
that  when  two  codewords  (corresponding  to  polynomials  of  degree  1)  of  the  basis  are  added,  two 
possibilities  exist:  Xi+Xj  ( j  >  i  + 1  mod  to)  andxi  +  Xj+i.  For  Xj+Xj,  {ei,e^_ i)(j), 

{ej,  ey-UG),  e(J')(j+1)}  =  {}  gives  the  location  of  overlapping  ones,  so  the  Hamming  weight  is 
3  +  3  =  6.  For  x*  +Xj+i,  {ej,  0  {ej_i,  e(j_2)(i_i),  !)(*)}  =  {e^-r)^)}  gives 

the  location  of  overlapping  ones,  so  the  Hamming  weight  is  3  +  3  —  (2)(1)  =  4.  Both  cases 
yield  a  Hamming  weight  greater  than  3.  Additionally,  if  3  or  more  codewords  (corresponding 
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to  polynomials  of  degree  1)  are  added  (xi  +  Xj  +  Xk  +  ■  •  ■) ,  then  e,,  e^,  e*, . . .  are  not  in  the 
intersection,  so  there  are  at  least  3  non-overlapping  ones. 

Now,  care  must  be  exercised  in  adding  the  codeword  (1,1,.. .,1).  When  this  codeword 
is  added  to  any  other  codeword  with  Hamming  weight  w,  the  Hamming  weight  becomes 
2 m  +  1  —  w.  Note  that  when  in  =  3,  when  two  codewords  corresponding  to  polynomials 
of  degree  1  are  added,  the  case  is  never  Xi  +  Xj ,  making  the  Hamming  weight  4.  Adding 
any  two  codewords  (corresponding  to  polynomials  of  degree  1)  with  the  constant  polynomial 
results  in  2m  +  1  —  4  =  7  —  4  =  3.  For  m  >  4,  2m  +  1  —  6  >  (2) (4)  —  5  =  3  (the  Xi+Xj  case) 
and  2 m  +  1  —  4  >  (2) (4)  —  4  =  5  (the  Xi  +  Xi+\  case),  so  the  minimum  Hamming  weight  is 
still  3  after  any  two  codewords  (corresponding  to  polynomials  of  degree  1)  are  added  with  the 
codeword  (1,1,...,  1).  When  3  or  more  codewords  (corresponding  to  polynomials  of  degree 
1)  are  added,  two  cases  arise:  either  3  or  more  points  evaluate  to  1  on  multiple  polynomials, 
or  fewer  than  3  points  have  this  overlap.  If  the  case  is  the  former,  then  adding  the  sum  with 
the  codeword  (1,1,...,  1)  yields  a  codeword  with  Hamming  weight  of  at  least  3  (since  the 
overlapped  points  had  ones  that  added  to  0  since  only  two  codewords  [x»,  Xi+ 1]  can  overlap  at 
e(*)(*+i))-  If  only  2  points  overlap,  then  it  is  the  case  Xi-i  +  Xi  +  Xi+ 1  +  . . . ,  with  Xi+2  ^  Xi-i 
(or  else  3  points  would  overlap).  Note  then  that  Xi+2  cannot  be  in  the  sum  of  codewords 
(since  it  would  evaluate  to  one  on  the  e(,;+i)p+2)  with  Xi+ 1,  creating  an  additional  overlapping 
point).  Thus,  the  component  of  the  sum  corresponding  to  e*+i  is  0.  When  the  codewords 
are  added,  the  components  corresponding  to  the  overlapping  points  become  zero  (since  only 
two  codewords  overlap  at  those  points),  in  addition  to  the  component  corresponding  to  ei+i . 
Thus,  when  (1,1,...,  1)  is  added,  at  least  3  components  of  the  vector  are  1,  so  the  Hamming 
weight  is  at  least  3.  In  the  case  with  1  point  overlapping,  the  sum  must  be  Xi  +  Xi+\  +Xj  + . . . 
with  Xi- 1  and  Xi+2  not  in  the  sum.  Thus,  the  components  corresponding  to  ej_i  and  e,;+i  are 
0,  along  with  (because  of  the  overlap).  Hence  when  the  sum  is  added  to  (1,1,...,  1), 

the  Hamming  weight  is  at  least  3.  If  no  points  are  overlapping,  then  the  sum  of  codewords  is 
Xi  +  Xj  +  Xk  +  ■  ■  ■  with  Xj+ 1,  Xj+ 1,  and  Xk+i  all  not  in  the  sum  (clearly,  Xi+i  ^  Xj+ \  ^  x/c+i). 
Thus,  in  the  sum  of  these  codewords,  the  components  corresponding  to  ej+i,  e7+i,  and  eu+i 
are  0.  Therefore,  when  the  sum  is  added  to  (1, 1, . . . ,  1),  the  Hamming  weight  is  at  least  3. 
Thus,  no  matter  how  many  codewords  are  added,  the  Hamming  weight  remains  at  least  3. 
Therefore,  d=  3.  □ 

Proposition  6.10.  The  dimension  of  C(Apm,  2)  is  2m  +  1,  and  the  minimum  distance  is  1. 

j  2 

Proof.  Recall  that  k  =  Yf  /,;.  Thus,  k  =  Y  fi  =  /o  +  fi  +  fi  —  1  +  m  +  m  =  2m  +  1.  By 
2  —  0  2  —  0 

Proposition  6.2,  n  =  2m  +  1.  Thus,  the  generating  matrix  for  the  code  can  be  made  into  an 
upper  triangular  (2?n  +  1)  x  (2m  +  1)  matrix  with  Is  along  the  diagonal  (since  the  row  vectors 
are  linearly  independent).  Thus  the  last  row  of  the  vector  must  have  have  all  nonzero  entries 
with  the  exception  of  the  last  entry  (that  is,  (0, 0, ... ,  0, 1)),  giving  it  a  Hamming  weight  of  1. 
Therefore,  the  minimum  Hamming  weight  (and  the  minimum  distance)  of  the  codes  is  1.  □ 

Remark  6.11.  The  general  characterization  of  binary  codes  based  upon  a  simplicial  complex 
which  is  an  m-gon  is  as  follows.  The  case  with  j  =  0  is  trivial  and  not  included.  Work  above 
indicates  that  if  j  =  1,  a  [2m  + 1,  m  + 1, 3]2  code  results,  whereas  if  j  >  2,  a  [2m+l,  2m+l,  1]2 
code  results.  Note  that  if  j  >  3,  no  more  linearly  independent  row  vectors  are  created  (since 
there  are  no  points  in  the  ideal  with  three  components  that  are  1  since  only  edges  are  vertices 
are  present  in  the  simplicial  complex). 

6.3.  Skeletal  Codes.  Computing  the  minimum  distance  for  a  general  simplicial  complex 
evaluation  code  at  this  time  seems  difficult.  In  this  section  we  focus  on  one  case  where  the 
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minimum  distance  can  be  calculated.  The  proof  is  elementary  but  slightly  involved.  To  begin 
we  define  the  codes  and  examine  a  few  examples. 

Definition  6.12.  For  0  <  h  <  £  an  h-skeleton  is  a  simplicial  complex,  denoted  by  A{£,  h),  on 
£  vertices  consisting  of  all  possible  h  —  1  to  O-dimensional  faces. 

Definition  6.13.  A  binary  h-skeleton  code  is  the  binary  evaluation  code  C(AAU,h)i3)i  °f  the 
associated  coordinate  arrangement  to  the  /i-skeleton  and  is  denoted  K(£.  h,j). 

The  length  of  these  codes  is  obtained  by  counting  the  number  of  points  that  have  at  most 
h  non-zero  entries.  Hence,  the  length  of  K(£,h,j)  is 


The  dimension  can  be  found  as  a  specialization  of  Proposition  6.2. 

Corollary  6.14.  For  0  <  j  <  h  <  £ 

dim  (K(£,h,j))  =  ^ 

i=0 

For  the  remainder  of  this  section  we  study  the  minimum  distance  of  the  codes  K(£ ,  h,j).  In 
order  to  obtain  more  information  about  K(£,  h,j)  we  need  to  examine  and  carefully  construct 
a  convenient  generating  matrix.  To  do  this  we  need  a  little  notation.  Let  a  =  {*i, . . . ,  ir}  C  [£\ 
and  let  xa  =  ■  ■  ■  Xir.  With  this  notation,  the  Stanley-Reisner  ideal  of  A(£,  h )  is 

I A  (l,h)  =  (Xa  ■  M  =h+  1). 

To  denote  points  in  the  arrangement  AA(e,h)  we  let  {e\ . . . ,  ei }  be  the  standard  basis  for  V  =  Fj 
(that  is,  e*  has  all  components  0  except  a  1  in  the  *-th  coordinate).  Now  for  r  =  { i\, . . . ,  *s}  C 

S 

[£\  let  eT  =  ^  eik.  Then  the  set  of  all  points  in  the  arrangement  AA(i,h)  is 
k= 1 


u  x 


=  {eT  :  1  <  |r|  <  h}. 


Now  we  will  construct  the  blocks  of  the  generating  matrix  for  I\(£,h,j).  Let  Brs  be  the 
matrix  defined  as 


(2)  Brs  =  (xa(eT)) 

where  \a\  =  r,  |r|  =  s,  and  the  rows  and  columns  are  ordered  lexicographically.  Then  a 
generating  matrix  G(£,h,j)  of  K(£,h,j )  constructed  block- wise  is 

G(£.  h,  j)  =  (Hrs)o<r<j- 

0<s</i 

We  can  now  denote  column  and  row  blocks  of  the  generating  matrix. 

Definition  6.15.  Let  CBt  =  {Brt  :  0  <  r  <  h}  be  the  union  of  the  blocks  of  columns  in  the 
matrix  of  the  code  with  t  ones  in  each  point.  Let  RBt  =  { Btr  :  0  <  r  <  j}  be  the  union  of  the 
blocks  of  rows  in  the  matrix  of  the  code  with  t  ones  in  each  point. 


This  notation  for  the  generating  matrix  streamlines  the  computation  of  minimum  distance. 
We  begin  by  presenting  an  upper  bound  for  the  minimum  distance. 
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Lemma  6.16.  For  0  <  j  <  h  <  £  the  minimum  distance  of  K(£,  h,j)  satisfies 


h-j 

E 


i- 0 


Proof.  In  the  generating  matrix  G(£,  h,  j)  the  rows  in  the  last  row  block  RBj  have  the  smallest 
weight.  The  smallest  t  such  that  Bjt  has  no  zero  entries  is  when  t  =  j.  The  Hamming  weight 
of  any  row  of  Bjt  for  j  <  t  <  h  is  (*j~ j) .  Hence,  the  Hamming  weight  of  an  entire  row  in  RBj 
is 


E 

i=j 


□ 


Remark  6.17.  If  j  =  h  then  K(£,  h,j)  is  a  maximum  distance  separable  (MDS)  code  but  the 
minimum  distance  is  1  because  the  upper  bound  here  is  1.  An  MDS  code  is  one  that  satisfies 
the  Singleton  bound  of  k  <  n  —  d  +  1. 

h-i 

If  j  =  1  then  the  minimum  distance  is  bounded  by  d  <  (E)-  ^  seems  difficult  to 

i=o  * 

determine  the  actual  minimum  distance.  The  main  result  of  this  paper  (Theorem  6.22)  is  that 
this  upper  bound  is  exactly  the  minimum  distance  for  the  case  j  =  1.  First,  we  obtain  a 
formula  for  the  Hamming  weight  of  adding  s  rows  of  the  generating  matrix.  In  order  to  do 
develop  this  formula,  we  need  a  little  more  notation.  Suppose  , . . . ,  Xis  are  the  degree  one 
monomials  that  correspond  to  the  s  rows  we  are  to  sum  in  B\a.  Let  Pa  be  the  set  of  all  points 
in  that  have  exactly  a  nonzero  entries.  Note  that  Pa  corresponds  to  the  columns  of  Hla. 

Definition  6.18.  For  1  <  t  <  s,  let  Xr  :=  {p  £  Pa  :  Xir(p)  =  1}.  Let  £“’s  be  the  set  of  all 
the  sets  of  points  that  evaluate  to  1  on  at  least  t  degree  one  monomials,  so 

cfi3  =  {Xkl  n  •  •  •  nllf  :  {fci, . . .  ,kt}  c  {*!, . . .  ,is}}. 

If  we  wanted  to  calculate  the  size  of  the  union  of  the  sets  X,,  U  •  •  •  U  Xis ,  then  we  could  use 
a  standard  inclusion-exclusion  formula 

S 

ix,;iu-..uxj  =  ^(— ir1  £  m 

t=i  re-cr 

However,  we  want  to  calculate  the  Hamming  weight  of  the  sum  of  these  row  vectors  of  which 
not  all  points  will  sum  to  1.  To  do  this  we  will  create  a  generalized  inclusion-exclusion  formula. 

Lemma  6.19.  The  Hamming  weight  of  adding  s  row  vectors  of  B\a  of  the  code  C(AA(i,h),  1) 
is 

E(-2 y-1  £  in 

i=i  Yecps 

Proof.  We  prove  this  by  induction.  The  critical  idea  here  is  that  if  a  point  p  is  contained  in 
exactly  t  sets  Xkl , . . . ,  Xkt  and  not  in  any  others,  then  the  entry  corresponding  to  this  point  in 
the  sum  will  be  0  if  t  is  even  and  1  if  t  is  odd.  Let  c*  be  the  coefficient  that  will  be  multiplied 
to  the  point  p  that  is  contained  in  exactly  t  sets  Xkl , . . . ,  Xkt  in  the  sum  (note  that  points  are 
the  objects  being  summed  here  because  the  Y s  consist  of  points).  In  the  case  when  t  =  1,  we 
want  to  count  all  the  points  that  are  in  exactly  1  set.  We  therefore  sum  the  entirety  of  the  sets 
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of  just  one  intersection:  ^  Xir  \ .  Thus,  the  coefficient  is  Ci  =  1  for  the  t  =  1  term.  However, 

r=  1 

if  t  >  1,  the  point  p  has  already  been  counted  in  lower  terms  because  it  is  also  a  subset  of  all 
possible  intersections  of  these  t  sets: 

xkl ,  •  ■  • ,  xkt ,  xkl  n  xk2 , . . . ,  xktl  n  xkt , . . . ,  xkl  n  •  •  •  n  xkt_1 , . . . ,  xk2  n  •  •  •  xkt . 

Because  we  want  the  coefficient  for  t  odd  to  be  1  and  for  t  even  to  be  zero,  we  now  have  that 

t 


E 


t,  odd 
t  even 


One  method  to  do  this  is  to  set 


E 

r—  1 


Cp  — 


(-!)*-! 
-2 


Now  we  prove  by  induction  on  t  that  Ct  =  (—2 )*  1 .  The  base  is  already  provided  above. 

/  *j.  \  t  “I- 1  ^  ^  ^  .  i  j  i . 

By  construction,  ct+i  =  _2 - (  r  )Cr"  Then  by  the  induction  hypothesis, 


Ct+l  — 


(-l)t+i  -  1  *{t  +  1 


-2 


-E 


(-2) 


r—  1 


Using  the  binomial  expansion  formula,  we  see  that  that 


t+i 


2=0 


(-l)t+1  =  (-2  +  l)t+1  =  E  ( 1  +  1 )  (-2)ilt+1_i  = 


t+1 


*=0 


t+1 


(—2)* 


t+ 1 


=  1  +  U('t1)(-2>‘  =  1  +  (-2>E 


t+1 


t+1 


(-2  y 


Hence, 


(-l)t+1-f 


-2 


=  E  7  (-2r1  =  (-2)t  +  E 


i= 1  x  7  i—1 

Now  add  the  sum  to  both  sides  of  this  equation  to  obtain 


t+1 


(—2  y 


(-2  Y  = 


(—!)*+!-!  *  /i  +  i 


—2 


-E 


(-2)i-1=ct. 


□ 


Lemma  6.19  gives  a  nice  method  to  compute  the  Hamming  weight  of  the  sum  of  s  row 
vectors. 


Lemma  6.20.  The  Hamming  weight  of  adding  s  vectors  in  RB\  is 

h  s 

£B-2>'-1 


e-t 


tj  \a  —  t 
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Proof.  Note  that  for  any  Y  £  Cf’s ,  the  size  is  |Y|  =  (^_*)  since  t  of  the  nonzero  entries  must 
match  up  with  the  t  monomials  (so  there  are  a  —  t  ones  left  to  chose  from  the  remaining  i  —  t 
components  of  the  point).  Since  we  can  choose  any  t  subsets  of  the  monomials,  we  have 


ici  = 


Then  the  formula  for  the  Hamming  weight  is  given  by  applying  Lemma  6.19  and  summing 
over  all  possible  column  blocks  CBa  where  1  <  a  <  h.  □ 

Next  we  prove  a  technical  lemma  that  will  be  used  in  the  proof  of  the  main  theorem. 
Lemma  6.21.  If  g?  =  J2  (®I*)  (®)  (-2)t_1,  then 


t=  1 


9i  = 


_  /  0  2|i 


0  2 \i  ' 


Proof.  We  prove  this  in  two  cases.  Case  1  is  when  i  =  2m.  Then 


9  2r. 


=  E 


/  s  —  t 
\2?n  —  tj\t 


=  V  (*-*)!*! 

^  (2m  —  t)\(s  —  2m)!f!(s  —  t)\ 


(-2) 


t- 1 


=  E 


si 


(2m)! 


— j'  (2m  —  t)!(s  —  2 m)\t\  (2m)! 

2  m 


(— 2)t_1  = 


5! 


(s  —  2m)!(2m)!  (2m  —  t)\t\ 


E 


(2m)i  (_2)t-i 


5 

2m 


E 


s 

2m 


t=i 
2m 


2m 

t 


(-2 )‘ 


5 

\2m 


*=i 
2m 


E 

Lf=0 


/2m 

V  t 


(~ 2)*  —  1 


E 

.t— o 


/  2m 

v  t 


(— 2)4(1) 


t  ( ~\  \2m—t 


- 1 


s 

2?n 


)(4)  I*1-2)2”-1] 


=  (s 
\2m 

Case  2  is  when  i  =  2m  +  1.  Then 

2m+l 

t=l 

2m+l 


-2)d-l)=0. 


s  —  t 

2  m.  +  1  —  tj  \  t 
( s  —  t)\s\ 


(2m  +  1  —  t)\(s  —  2m  —  l)!f!(s  —  £)! 
2y^1  s!  (2m +1)! 


(— 2)t_1 
(-2)4"1 
(-2)4-1 


— (  (2 m  +  1  —  f)!(s  —  2m  —  l)!t!  (2m  +  1)! 

V1  (2m +  1)! 


(s  —  2m  —  l)!(2m  +  1)!  (  (2m+l  — t)!t! 

f2m+l 


S 

2?n  +  1 


s 

2m  +  1 


E 


t= l 
'2m+l 

E 

.  t=0 


2m  +  1 
t 


2m  + 1 
t 


{-2f 

(-2)4-l 
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2  m  +  1 


'2m +1 

E 

.  t= o 


/  2  m  +  1 

v  t 


(— 2)*(1) 


2m+l— t 


- 1 


s 

2?n  +  1 


[(1-2)"-!]  = 


5 

2m 


-5  )  (-2>  = 


f s 

\2m 

(’ 

\2m 


-2K-1-1) 


□ 


Notice  that  the  formula  in  Lemma  6.20  for  the  case  s  =  1  is  exactly  the  computation  made 
in  Lemma  6.16.  In  order  to  show  that  this  value  for  s  =  1  is  exactly  the  minimum  distance,  it 
is  enough  to  show  that  the  sum  for  s  >  1  is  greater  than  that  for  s  =  1,  since  q  =  2.  Now  we 
can  state  and  prove  the  main  theorem  of  this  paper. 

Theorem  6.22.  The  minimum  distance  of  the  codes  C(A/±uth)A)  —  K(t,h,  1)  is 

£(::!> 

a— 1  x  7 

Proof.  We  need  to  show  that  the  Hamming  weight  of  adding  s  rows  given  in  Lemma  6.20  is 
always  larger  than  the  Hamming  weight  of  one  row  given  in  Lemma  6.16: 


EE(-2>‘_‘ 

a=  1  t—1 

This  is  equivalent  to  showing 
h 


i-t 

a  —  t 


>E 


- 1 

a  —  1 


(3) 


E 

a— 1 


E'-2> 


t-i 


vt= 1 


£-  1 

a  —  1 


>  0. 


Now  we  use  Pascal’s  formula  to  allow  for  the  exchange  of  terms  of  3.  We  examine  the  term 


t-i 

a  —  t. 


’-t-  1 
a  —  t 


-t-  1 
a  — t—1 


'  -t-  2 
a  —  t 


t-t-2 
a  —  t—1 


t-t-2 

a  —  t—1 


—  t  —  3 
a  —  t  —  2 


Since,  whenever  Pascal’s  formula  is  used,  each  binomial  coefficient  is  broken  down  into  two 
binomial  coefficients,  the  process  is  analogous  to  Pascal’s  triangle:  the  top  number,  l  —  t  —  x, 
corresponds  to  the  xth  row,  and  the  bottom  number,  a  —  t  —  x,  corresponds  to  the  .rth  column. 

S 

Thus,  there  are  *)  occurrences  of  each  for  each  t.  Therefore,  *)  (i-t)  > 

so  the  inequality  we  are  trying  to  prove  is  now 
h 


i—  1 


(<)  E  E(-2)"' 

a= 1  L  \i=l  i= 1 

Now  focusing  on  the  ( tafs )  terms,  4  becomes 
h 

(5)  E 


>  0. 


E  Et-2)4'1 


.i= 1 


>  0 
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Notice  that  the  third  sum  is  only  nonzero  when  t  <  i  and  that  the  (s_j )  term  only  affects  the 
t  =  1  term  of  the  third  sum.  Hence,  we  can  rewrite  5  as 


t  +B-2)- 


(6)  E 

a=  1 

Let  df  be  the  coefficient  of  (*“®)  in  6: 


s\  I  s  —  t 


£  —  s 
a  —  i 


>  0 


Recall  the  numbers  gf  from  Lemma  6.21: 

9t=iz 


s  —  t\  (  s 


i  —  t)  \t 


Then 


s  —  ( s 


i  —  t]  \t 


(-2) 


s\  /  s  —  t 


t- i 


s \  ( s  —  t 


t  J  \i  —  t 


s  —  1 
i  —  1 


Thus, 

(7)  ^2m+l  =  92m+ 1 

which,  by  Lemma  6.21,  gives  that  7  becomes 

d: 

Using  Lemma  7  again,  we  get 
4m  =  92m  ~ 


S-i\ 
2m  ' 


x2m+l 


5  —  1 
2m 


( s 

\2m 


s  —  1 
2  to  —  1 


s  —  1 
2m 


s  —  1 
k2  m  —  1 

Hence, 

(8) 

The  main  inequality  we  are  trying  to  prove,  6,  can  now  be  written  as 

h 


rls  —  —  ds 
u2m+l  “  u2m* 


(9) 


E 


E^s 

,i=i 


t  —  s 
a  —  i 


>  0. 


Assume  s  =  2m  is  even  and  expand  the  left  hand  side  of  9  via  odds  and  evens: 


E 


Ed; 


'2m 
2  r 


£  —  2m 
a  —  2  r 


(  m—  1 

Erfi 

V  r= 0 


2m 

2r+l 


Then  using  9  on  the  odd  terms,  we  get 
h 


(10) 


E 

a— 1 


Erfi 

Vr=l 


2m 

2r 


^  —  2?n 
a  —  2r 


-  E* 


'2  m 
2r 


s,  r— 0 


£  —  2m 
a  —  2r  —  1 


£  —  2m. 
a  —  2) —  1 
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Figure  9.  A(5, 2) 


Then  using  Pascal’s  formula  on  10,  we  have 


Then  switch  sums  on  11  to  get 


Then  the  sum  in  the  left  portion  of  12  telescopes: 


Then  notice  that  d^m  =  0  and  that  13  becomes 


Since  14  is  the  left  hand  side  of  6  and  each  term  is  positive,  we  have  proved  the  theorem.  □ 


Example  6.23.  A(5,  2)consists  of  5  vertices  and  all  1-dimensional  and  0-dimensional  faces 
(which  are  the  edges  and  vertices,  respectively).  This  object  is  a  4-dimensional  object,  but  it 
can  be  approximated  in  Figure  9.  The  matrix  generating  K (5,  2, 1)  is  as  follows: 

'1111111111111111" 
0100001111000000 
0010001000111000 
0001000100100110  ■ 
0000100010010101 
0000010001001011 


Thus,  by  the  formulas  given  in  Equation  1,  Corollary  6.14,  and  Theorem  6.22,  K{ 5,2, 1)  is  a 
IE  (').  E  (').  E  (1:1)]2  =  [E  (“)■  E  (“)■  E  (":i)l2  =  [16,6.5],  code. 

2=0  2=0  a—1  i— 0  2—0  a— 1 
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Now  we  focus  on  the  case  when  j  >  1.  It  is  more  complicated,  and  we  are  not  able 
to  calculate  the  minimum  distance.  However,  we  are  able  to  find  formulas  for  summing  row 
vectors  of  the  generating  matrix  and  and  are  able  to  compare  these  formulas  to  the  conjectured 
upper  bound. 

Definition  6.24.  Let  B(a 1; . . . ,  crs)  be  the  Hamming  weight  of  adding  s  rows  of  the  generating 
matrix  G(i,  h,  j)  where  each  row  corresponds  to  a  set  cq  C  [£]  =  {1, . . .  ,£}. 

The  j  >  1  analogue  to  Definition  6.18  is  the  following. 

Definition  6.25.  For  1  <  t  <  s  let  X(Tr  :=  {p  G  Pa  :  x„r{p)  =  1}.  Let  be  the  set  of  all 
the  sets  of  points  that  evaluate  to  1  on  at  least  t  of  the  s  monomials  xai , . . . ,  x„s .  Thus, 

Cat'8  =  {Xkl  :  {ki,...,kt}  c  {err, ... ,  crs}}. 

Proposition  6.26.  For  any  skeletal  code  K(l,h,j),  the  Hamming  weight  of  adding  n  rows  of 
the  generating  matrix  is 


n 

B(a1,...,crn)  =  ^2 

(— 2)e_1  53  B  | 

(l>) 

n 

=  E 

(-^•eT'C-T1) 

e=l 

IC[n\ 

L  \I \=e 

Kiel  J 

e=l 

I C  [n]  i= 0  V  J 

L  m=e  J 

Proof.  The  (— 2)t_1  coefficient  follows  with  the  same  argument  as  in  Lemma  6.19.  Then  we 
realize  that  the  number  of  points  in  a  f-fold  intersection  Y  =  Xkl  D  •  •  •  D  Xkl  £  Cf’s  is  equal 
to  the  Hamming  weight  of  the  row  corresponding  to  the  union  of  the  sets 

t 

2—1 

Note  that  this  row  might  not  actually  exist  in  the  generating  matrix.  However,  we  can  consider 
it  as  a  row  in  the  full  matrix  where  j  =  i.  Finally,  the  remainder  of  the  formula  is  realized  by 
applying  Lemma  6.16.  □ 


In  order  to  show  that  the  minimum  distance  is  equal  to  the  upper  bound,  it  must  be 
demonstrated  that 

h—  |Uf7i  | 

<-2r‘E  E 

/C[n]  2=0 

m=e 

However,  this  proof  turns  out  to  be  difficult.  For  the  remainder,  we  examine  a  few  cases. 

h—j 

Proposition  6.27.  For  j  >  1,  .B(<7i,  (j2)  >  E  (V)' 

2=1 

Proof.  Let  |<7i|  =  i\,  |cx2 1  —  ii-  Without  the  loss  of  generality,  assume  (relabeling  as  necessary) 
i\  <  f 2  and  (T i  ^  o 2  (else,  B{a\1a2)  =  0).  Note  that  |cri  U  CT2 1  >  i\  +  1.  By  Proposition  6.26, 

=  B(a1)  +  B(a2)  -  2B(cti  U  a2) 


(15) 


E 

i—1 


e-j 

i 


^E 


h—ii 

=  E 


i=  1 


-  h 
i 


h—i  2 

+  E 


i= 1 


h—  |<tiU<T2| 

2  E 

i—1 


(e  -  | (7\  U  cr2|^ 
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Since  |  cri  U  cr2 1  >  ii  +  1, 


h—ii 

E 


2= 1 


-  h 
i 


h—i  2 

+  E 


2= 1 


/l— |(TiU(J2| 

2  E 

2=1 


^  -  |<Ji  U  cr2|^ 


h— i\ 

>E 


h—ii 

+  E 


-  ^2 
i 


h—  (ii  +  1) 

2  E 


(ll  + 1) 


2—1  N  '  2=1  N  '  2=1 

Hence,  we  prove  the  inequality  for  \<7\  U  <r2|  =  *i  +  1.  There  are  now  two  cases:  (1)  i\  =  z2 
or  (2)  i\  +  1  =  z2  and  o\  C  cr2.  For  Case  1,  assume  Zi  =  z2  =  a,  so  the  left-hand  side  of  the 
proposed  inequality  becomes 


Since  j  >  a, 


(E)  —  E  C  i°)i  so  it  suffices  to  show 

i- 1  i=  1 


This  last  expression  can  be  rewritten  as 


We  can  phrase  this  last  inequality  by  saying,  “The  Hamming  Weight  decreases  by  more  than 
half  in  each  group  of  rows  (each  group  of  rows  corresponds  to  polynomials  of  the  same  degree).” 

Since  ©  =  (Y)  +  (rl). 


h—j 

Thus,  for  Case  1,  B(cti,ct2)  >  J2  (E)‘ 

2=1 

For  Case  2,  since  the  Hamming  Weight  decreases  by  more  than  half  in  each  group  of  rows, 
B{a1}  a2)  >  B(a1)  -  B(a2)  >  E(cti)  >  B(a2)  >  B(aj), 


where  B{(jj) 


E  (E)>  since  \ai\  -  l°"2|-  Thus,  B{al,a2)  >  £  (E)‘ 


□ 
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Proposition  6.28.  B(ai, . . . ,  an)  =  B(a1, ...,  +  B(an)  -  2B(ai  U  u„, . . . ,  ct„_i  U  <t„). 

Proof.  From  Proposition  6.26, 

n 

(16)  B(ai,  ...,an)  =  'Yl  B(cri )  -  2  ^  B(aii  u  CTi2)  d - P  (-2)n_1B(cri  U  •  •  •  U  ern). 

i— 1  ii,22 

Rearranging  so  that  all  terms  involving  crn  on  the  right-hand  side  of  16  are  isolated  yields 


n—  1 

(17)  5Z  %U<7i2)  +  ...  +  (-2)"-2%U-Uff„_i) 

£=1  H  >£27^™ 

n— 1 

B(o’n)  —  2  ^  '  73(ci  U  <t„)  +  •  •  •  +  (—2)"  1B((Ji  U  •  •  •  U  crn) 

i=l 

Again,  by  Proposition  6.26, 

n—  1 

(18)  B(o-i,...,qn_i)  =  y^E(o-j)-2  ^  B(aix  U  ai2)  -f - h  {-2)n~2B(ai  U  •  •  •  U  cr„_i), 

i= 1  £1,127^ 

and,  since  cr,  U  cr,  U  cr„  =  (cr*  U  <t„)  U  (crj  U  <jn), 

(19)  B(cti  Ucr„, . . .  ,<T„_ 1  U  <rn) 

n—  1 

=  J2  B((Ji  U  cr„)  -  2  )T)  B(ail  U  cri2  U  an)  H - b  {-2)n~2 B{ai  U  •  •  •  U  <r„). 

i=l  »i,i2 

Substituting  18  and  19  back  into  the  right-hand  side  of  17  yields  the  claimed  formula.  □ 

Proposition  6.29.  Let  |oi|  <  |cr2 1  <  •  •  •  <  |crfc|  <  |crfe_|_i |  =  •  •  •  =  |cr„_i|  =  \an\  for  k  <  n.  If 
1 0*1  <  (n  -  k)  +  |tr„|,  then  Bf Ji,  ...,an)>  B(an). 

Proof.  The  proof  will  be  by  induction  on  n.  The  base  case  n  =  2  is  covered  by  the  proof  of 
Proposition  6.27.  Thus,  assume  B(a±, . . . ,  ax)  >  B(ax)  for  x  <  n.  It  must  be  demonstrated 
that  B(<ji,  . . .  ,an)  >  B(an).  By  the  formula  in  Proposition  6.28, 

B(<Ji,.  ■ .  ,<rn)  =  B(ai, . . .  ,0n- 1)  +  B(an)  -  2B(ai  U  <j„,  . . .  ,ct„-i  U  ct„). 

Note  that  B(ai  U  an, . . . ,  cr„_i  U  an)  is  the  Hamming  weight  of  n  —  1  row  vectors  being  added 
together,  fulfilling  the  property  that  all  nonzero  entries  of  the  rows  Cj  U  an  are  also  nonzero 
entries  in  the  row  an.  Hence, 

B(an)  >  B(ai  U  an, . . . ,  ct„_  1  U  an). 

Thus, 

B(ai, . . . ,  cr„)  >  H(cti,  . . . ,  an- 1)  +  B(crn)  -  2 B(an) 

=  H(cti,  . . .  ,cr„_i)  -  B(an). 

Then  using  Proposition  6.28  repeatedly,  we  have 

B{( Ji, . . .  ,cr„)  >  B(<ji,  . . .  ,cr„_2)  +  B(an- 1)  -  2B(ai  U  er„_i, . . . ,  cr„_2  U  cr„_i)  -  B(an) 


>  B(a  1, . . .  ,cr„_2)  +  B(an- 1)  -  2B{an-1)  -  B{an) 

=  B(a1, . . . ,  cr„_2)  -  5(crn_i)  -  B(an) 

>  ■■■>  B(ai, . . . ,  crfc)  -  (B(<jk+i)  H - h  B(crn)). 
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By  the  induction  hypothesis, 

B(ai, . . .  ,an)  >  B(ai , . .  ,,ak)  -  (B(crk+ 1)  H - b  B(crn)) 

>  B(ak )  —  (B(crk+i)  +  •  •  •  +  B(an)) 

=  B(ak)  -  (n-  k)B(an), 

since  |erfc+i|  =  •  •  •  =  |er„|.  Let  Tk+m  C  {£]  such  that  |Tfc+m|  =  |crfc|  +  m  (so  B(a„)  <  B(Tk+m) 
for  to  <  (n  —  fc))  and  recall  from  Proposition  6.27  that  the  Hamming  weight  between  groups  of 
rows  decreases  by  more  than  half  (also,  note  that  Tk+i  is  a  row  with  higher  Hamming  weight 
than  the  row  Tk+i+i).  Thus, 

. . . ,  a„)  >  B(ak)  -  (n  -  k)B(crn) 

>  [ B{ak )  -  B(Tk+ 1)]  -  B(Tk+2 ) - B{an) 

>  [B{Tk+i)  -  B(Tk+2 )]  -  B(Tk+3) - B(an) 

>  . . .  >  B{Tn_i)  -  B(<jn )  >  B(an). 

Note  that  there  are  exactly  enough  Tk+m  terms  for  the  above  inequalities  since  |  cr^  |  <  (n  — 
k)  +  \an\.  Thus,  for  |crfc|  <  (n  -  k)  +  \an\,  B(ai, . . .  ,tr„)  >  B(an).  □ 


With  respect  to  Equation  15  Propositions  6.26,  6.27,  6.28,  and  6.29  give  evidence  for  the 
following  conjecture. 


Conjecture  6.30.  The  minimum  distance  of  K(£,h,j)  is 

h-j 

d  =  E 


i- 1 


J 


Proposition  6.31.  The  codes  K(£,£—  l,t—  2)  and  He  are  permutation  equivalent. 

Proof.  By  Theorem  1.8.1  in  Huffman  and  Pless  [9],  any  [2r  —  1,  2r  —  1  —  r,  3]  code  is  equivalent 
to  Hr-  The  length  of  K (£,  l  —  1,  £  —  2)  by  Equation  1  is 


i= o 


e-i 


i=0 


£  =£  -E  - 


i= 0 


=  2^  —  1. 


The  dimension  of  i  —  1,^  —  2)  by  Corollary  6.14  is 


t- 2 


1—2 


e/.=e/.=e  :  -e  - 


0  i= 0  i=0  x  7  2—0 

The  upper  bound  on  the  minimum  distance  Lemma  6.16  is  by 


h-j 

d<E 


j 


e-(e-i) 

=  E 


i= 0  x  '  y  i= 0  x  y  i=0 


£  -  (£  -  2) 


=  E 


=  3. 


By  Corollary  1.4.14  in  [9],  since  the  rows  in  the  matrix  generating  K(£,£—  l,£  —  2)  are  linearly 
independent,  d  >  3,  implying  d  =  3.  Thus,  I\(£,£  —  1,£  —  2)  =  He-  □ 


Remark  6.32.  Since  the  two  codes  are  permutation-equivalent,  there  must  be  some  means 
by  which  to  transform  K(£,£  —  1,£  —  2)  into  He-  This  process  is  achieved  by  adding  rows, 
permuting  columns,  and  permuting  rows  in  the  matrix  generating  K (£,  £  —  1,  £  —  2).  Since  the 
matrix  (if  the  points  in  V (/)  are  ordered  appropriately)  can  be  made  into  an  upper-triangular 
matrix  adjoined  to  a  (2f  —  £  —  1)  x  £  matrix,  adding  all  rows  to  rows  above  them  will  result  in  a 
[I2£ -g-i  1^4]  matrix.  In  the  A  block  (corresponding  to  points  in  V{I)  that  have  only  one  0  [and 


SUBSPACE  ARRANGEMENT  CODES  AND  CRYPTOSYSTEMS 


39 


thus  l—l  ones]  as  a  component),  there  are  (since  there  are  £—2  rows)  ]T)  (  f  3)  =  —  1 

ones,  which  is  odd.  Thus,  when  all  the  rows  are  added,  the  A  block  remains  unchanged.  This 
A  block  consists  of  binary  strings  of  length  £  that  corresponding  to  polynomials  up  to  degree 
£  —  2  being  evaluated  on  points  with  £  —  1  ones.  These  strings,  once  transposed  and  combined 
with  a  transposed  identity  matrix,  will  consist  of  all  numbers  in  binary  from  1  to  2e,  indicating 
its  equivalence  to  He. 


7.  McEliece  Cryptosystem 

An  important  application  of  the  working  on  coding  theory  and  subspace  arrangement  codes 
is  the  creation  of  a  cryptosystem.  Work  has  already  been  done  on  making  ciphers  based  upon 
error-correcting  codes.  The  most  prominent  example  is  the  McEliece  Cryptosystem,  which  is 
described  in  detail  in  [6]  and  [18]. 

7.1.  History  and  Description.  Robert  McEliece  first  proposed  the  idea  of  a  public  key 
cryptosystem  built  from  error-correcting  codes  in  1978.  Since  McEliece’s  development  of  the 
scheme,  various  proposals  have  been  made  with  different  error-correcting  codes,  but  most  of 
them  have  been  proven  to  not  be  as  efficient  as  McEliece’s  proposal  to  use  Goppa  Codes  [18]. 
Goppa  codes  have  a  closed  formula  for  length,  dimension,  and  minimum  distance  based  upon 
two  parameters,  m  and  t:  n  =  2m,  k  =  n  —  mt ,  and  d  =  2f  +  1.  The  strength  of  the  McEliece 
cryptosystem  is  that  it  can  be  very  hard  to  determine  what  codeword  is  closest  to  a  random 
string  of  digits.  The  process  for  encryption  and  decryption  follows  a  few  basic  guidelines.  The 
receiver  will  make  a  public  key,  and  the  sender  will  encrypt  their  message  with  the  key  so  that 
the  receiver  can  decrypt  it  (note  that  in  the  standard  cryptologic  terminology,  the  receiver  and 
sender  are  Bob  and  Alice,  respectively).  The  receiver  selects  the  generating  matrix  G  for  C,  a 
[n,k,d\  2  code,  so  G  is  a  k  x  n  matrix.  The  receiver  also  selects  S,  an  invertible  k  x  k  matrix 
over  F2,  and  P,  an  n  x  n  permutation  matrix.  The  receiver  then  calculates  G\  =  SGP.  While 
keeping  S,  G,  and  P  private,  the  receiver  publishes  G\  as  a  public  key.  The  sender  who  wants 
to  send  a  message,  x,  randomly  selects  a  vector,  e,  of  length  n  with  weight  t.  The  sender  then 
forms  y  =  xG\  +  e  and  sends  y  to  the  receiver.  The  receiver  then  needs  to  decrypt  y  to  find 
x.  First,  the  receiver  finds  y\  by  calculating  yP~x .  Using  the  error  decoder  for  C  (such  as 
a  parity  check  matrix)  the  receiver  finds  x±,  the  codeword  in  C  that  the  is  the  closest  to  y-\ . 
Next,  the  receiver  finds  Xq,  where  x^G  =  X\.  Finally,  the  receiver  recovers  x  by  calculating 
XoS1-1-  This  process  can  best  by  explained  by  an  example. 

Example  7.1.  Let  C  be  the  [15,  5,  7] 2  skeletal  code  corresponding  to  A"(4,  2,1).  The  matrix 
for  C  is 

'111111111111111' 
010001110001110 
001001001101101  . 
000100101011011 
000010010110111 
However,  to  make  calculating  a  parity  check  matrix  easier,  replace  the  first  row  with  the  sum 


of  all  five  rows: 

'  1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

1 

0 

G  = 

0 

0 

1 

0 

0 

1 

0 

0 

1 

1 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

1 

0 

1 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

1 

1 
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'1  0  0  0  0  ' 

110  0  0 

Let  S'  =  0  110  0, 

0  0  110 
0  0  0  1  1 

[010000000000000' 
001000000000000 
000100000000000 
000010000000000 
000001000000000 
000000100000000 
000000010000000 
and  let  P  =  000000001000000  . 

000000000100000 
000000000010000 
000000000001000 
000000000000100 
000000000000010 
000000000000001 
_100000000000000_ 
[010000111111000' 
011000000111111 
Thus,  Gi  =  SGP  =  101100011110001  . 

000110110011011 
000011011110110 
Now,  the  sender  wants  to  send  the  message 

x  =  [  1  1  0  1  0  ], 

with  e  randomly  chosen  as 

e  =  [  1  1  1000000000000]. 

The  sender  then  computes  y  =  xG\  +  e,  resulting  in 

y=  [110110001011100]. 

The  receiver  then  starts  decrypting  y ,  first  calculating  yi  —  yP~l  =  yPlA 

ooooooooooooooi 
100000000000000 
010000000000000 
001000000000000 
000100000000000 
000010000000000 
000001000000000 
000000100000000 
000000010000000 
000000001000000 
000000000100000 
000000000010000 
000000000001000 
000000000000100 
000000000000010 

=  [101100010111001]. 


=  1  101  1000101  1  100' 
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The  next  step  is  to  use  the  parity  check  matrix,  H ,  of  C  to  correct  y\  for  errors.  Since 
G  =  [4 1<3],  by  principle,  H  =  [-QT\In_e].  Thus, 


H  = 


1  1  1 
1  1  0 
1  1  0 
1  0  1 
1  0  1 
10  0 
0  1  1 
0  1  1 
0  1  0 
0  0  1 


0  0  1 
1  0  0 
0  10 
1  0  0 
0  10 
1  1  0 
1  0  0 
0  10 
1  1  0 
1  1  0 


0  0  0 
1  0  0 
0  10 
0  0  1 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
1  0  0 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
1  0  0 
0  1  0 
0  0  1 


In  order  to  figure  out  how  to  correct  y\ ,  the  syndrome  of  y\ .  termed  S(yi),  needs  to  be  found. 

S(vi)  =  yiHT: 


101100010111001 


1  1  1 
1  1  1 
10  0 
0  1  0 
0  0  1 
10  0 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


1  1  1  0  0  0  0 
0  0  0  1  1  1  0 
110  110  1 
10  110  11 
0  110  111 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
1  0  0  0  0  0  0 
0  1  0  0  0  0  0 
0  0  1  0  0  0  0 
0  0  0  1  0  0  0 
0  0  0  0  1  0  0 
0  0  0  0  0  1  0 
0  0  0  0  0  0  1 


0  0  0  1  1  1  1  1  1  1 
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Note  that  the  syndrome  has  dimension  10,  so  there  are  210  =  1024  potential  error-correction 
actions  to  be  taken  (with  each  action  being  uniquely  correlated  to  a  syndrome).  Observe  that 

'1  1  1  1  1  10000' 

1  1  1  0  0  0  1  1  1  0 
10  0  110  110  1 
0  10  10  110  11 
0  0  10  110  111 
1000000000 
0100000000 

[110000000000001]  0010000000 

0001000000 

0000100000 

0000010000 

0000001000 

0000000100 

0000000010 

0000000001 

=  [0001  1  1  1  1  1  1], 

implying  that  [l  1  000000000000  l]  (which  is  equal  to  eP~l) 

is  the  coset  leader  for  the  syndrome  [000  1  1  1  1  1  1  1].  Thus,  to  error- 

correct  yi  to  find  X\ ,  add  the  coset  leader  to  y\ : 

xx  =  [101100010111001] 

+  [1  10000000000001] 

=  [011100010111000]. 

Next,  the  vector  XqG  =  X\  must  be  found.  Careful  observation  reveals  that  x\  is  the  sum 
of  the  second,  third,  and  fourth  rows  of  G,  so  Xq=  [  0  1  1  1  0  ] .  The  last  step  in  the 

process  is  to  compute  x  =  XoS1-1: 

'  1  0  0  0  0  ' 

110  0  0 

[01110]  11100  =[11010], 

11110 

11111 

The  original  message  has  been  correctly  deciphered. 

7.2.  Security.  An  important  reason  to  study  the  McEliece  cryptosystem  is  for  its  security. 
Work  done  by  Hang  Dinh,  Cristopher  Moore,  and  Alexander  Russell  in  [5]  indicate  that  Shor’s 
algorithm  is  not  effective  against  the  McEliece  cryptosystem.  Shor’s  algorithm,  discovered  in 
1994,  describes  a  factoring  process  possible  with  quantum  computing.  This  algorithm  can  be 
applied  to  modern  ciphers  as  a  quantum  Fourier  attack.  This  attack  has  been  demonstrated  to 
be  successful  on  a  number  of  cryptosystems,  including  the  widely-used  RSA  cipher.  However, 
Hang  Dinh,  Cristopher  Moore,  and  Alexander  Russell  show  in  [5]  that  Shor’s  algorithm  of 
factoring  integers  does  not  aid  in  solving  the  nearest  codeword  problem,  which  is  the  central 
point  of  security  of  the  McEliece  cryptosystem.  Hence,  the  McEliece  cryptosystem  is  immune 
to  this  attack.  Thus,  applying  skeletal  codes  to  the  McEliece  cryptosystem  could  result  in  a 
viable  method  for  encryption  if  quantum  computing  is  realized. 
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8.  Sage  Code 

An  integral  part  of  this  project  has  been  using  Sage  to  execute  calculations  that  would  be 
too  difficult  to  calculate  by  hand.  What  follows  are  four  script  sessions  that  were  part  of  the 
project. 

8.1.  Generating  Evaluation  Codes  and  Their  Properties.  Work  was  initially  done  in  the 
algebraic  geometric  software  Macaulay2.  The  following  code  (in  conjunction  with  an  example) 
details  using  a  Macaulay2  interface  in  Sage  to  construct  an  evaluation  code.  The  ideal,  finite 
field,  and  vector  space  are  used  to  create  the  matrix  representing  the  codes.  From  this  matrix, 
the  dimension  of  the  code  can  be  calculated. 


sage:  loadPackage  "RationalPoints" ; 

sage:  S=ZZ/7 [x,y ,z] ; 

sage:  I=ideal(x*y,x*z,y*z) ; 

sage:  p=rationalPoints  I; 

o3  :  Ideal  of 

sage:  codelength=length  p 

19 


sage:  p 

{-CO,  0,  0},  {1,  0,  0},  12 ,  0,  0},  {3,  0,  0},  {4,  0,  0>,  {5,  0,  0>,  {6,  0, 
0},  -CO,  1,  0},  -CO,  2,  0},  {0,  3,  Oh  -CO,  4,  0>,  {0,  5,  0},  {0,  6,  0>,  {0, 
0,  1},  {0,  0,  2},  {0,  0,  3>,  {0,  0,  4},  -CO,  0,  5>,  {0,  0,  6» 


List 

sage:  k=basis(0,4,S) 

I  1  x  x2  x3  x4  x3y  x3z  x2y  x2y2  x2yz  x2z  x2z2  xy  xy2  xy3  xy2z  xyz  xyz2  xz 


xz2  xz3  y  y2  y3  y4  y3z  y2z  y2z2  yz  yz2  yz3  z  z2  z3  z4  I 

1  35 

Matrix  S  < - 

sage:  e=apply(0  ..  (rank  source  k)-l,  i->(g_i=(a,b,c)->sub(k_(0,i) , 
{x=>a,y=>b,z=>c}))) ; 

sage:  T=matrix  toList  apply(0  ..  (rank  source  k)-l,j-> 
toList  apply (0  ..  #p-l,  i->e_j (toSequence  p_i) ) ) ; 

ZZ  35  ZZ  19 

o43  :  Matrix  ( — )  < - ( — ) 

7  7 


sage:  C=image  transpose  T; 
sage:  C 


image  1 1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

11 

1 

1 

1 

1 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

11 

2 

-3 

1 

2 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

11 

3 

2 

-1 

-3 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

11 

-3 

2 

1 

-3 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 

11 

-2 

-3 

-1 

2 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0 
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1 

-1 

1 

-1 

1 

0000000000000000 

0 

0 

0 

0 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

1 

1 

1 

1 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

2 

-3 

1 

2 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

3 

2 

-1 

-3 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

-3 

2 

1 

-3 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

-2 

-3 

-1 

2 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

-1 

1 

-1 

1 

000000 

0 

0 

0 

0  1 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

1 

1 

1 

1  1 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

2 

-3 

1 

2  1 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

3 

2 

-1 

-31 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

-3 

2 

1 

-31 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

-2 

-3 

-1 

2  1 

1 

0 

0 

0 

0 

0000000000000000 

0 

0 

0 

0 

000000 

-1 

1 

-1 

1  1 

ZZ  ZZ  19 

- module,  submodule  of  ( — ) 

7  7 

sage:  codedimension=rank  C 
13 


The  above  example  describes  an  evaluation  code  with  length  19  and  dimension  13  over  the 
finite  field  F7.  Note  that  Macaulay2  lacks  any  function  to  calculate  minimum  distance. 

The  following  code  uses  a  hie  assembled  through  help  from  Professor  W.  David  Joyner  that 
has  been  uploaded  to  the  Sage  server.  The  following  production  of  the  code  shows  the  various 
components  of  the  code  simultaneously  with  an  example. 

sage:  attach  /home/wakef ield/Downloads/evaluation_codes . sage 

sage :  R. <x,y ,z ,w>=PolynomialRing(GF (2) ,4,"x,y,z,w") ; 

sage:  I=R. ideal ( [x*y*z*w] ) ; 

sage :  P=get_exps (4,2) ; 

sage:  P 


[[0, 

0,  0, 

0], 

,  [1,  0,  0, 

0], 

,  [0, 

1, 

0, 

0], 

[0, 

0, 

1, 

0] , 

,  [0, 

0, 

0, 

1]] 

sage 

:  Pts 

=  get_pts (I) 

sage 

:  Pts 

[[0, 

0,  0, 

0], 

,  [1,  0,  0, 

0], 

,  [0, 

1, 

0, 

0], 

[1, 

1, 

0, 

0] , 

,  [0, 

0, 

1. 

0], 

[1, 

0,  1, 

0], 

[0,  1,  1, 

0], 

[1, 

1, 

1, 

0], 

[0, 

0, 

0, 

1], 

[1, 

0, 

0, 

1], 

[0, 

1,  0, 

1], 

[i,  1,  0, 

1], 

[0, 

0, 

1, 

1], 

[1, 

0, 

1, 

1], 

[0, 

1, 

1, 

1]] 

sage:  C  =  evaluation_code(P,  Pts) ;  C 

Linear  code  of  length  15,  dimension  5  over  Finite  Field  of  size  2 
sage:  C.minimum_distance() 

7 


The  first  line  of  code  attaches  the  hie  Professor  Joyner  created.  The  second  establishes 
the  hnite  held,  and  the  third  creates  the  ideal.  The  next  two  lines  find  the  points  on  the 
ideal  based  on  the  selection  of  j,  with  the  following  line  formally  establishing  the  code.  A 
built-in  Sage  function  for  minimum  distance  is  used  in  the  last  line.  The  example  uses  the 
ideal  <  X1X2X3X4,  >  in  F2  with  j  =  1  to  hnd  a  [15,5,  7]2  code.  An  important  note  is  that 
with  moderately-high  parameters  used,  this  Sage  code  is  not  efficient  and  can  take  signihcant 
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amounts  of  time,  especially  to  calculate  minimum  distance.  Thus,  it  is  very  beneficial  to 
develop  formulas  to  calculate  the  properties. 


8.2.  Boolean  Arrangement  Codes.  A  key  part  of  the  arguments  involved  with  developing 
formulas  for  the  dimension  of  error-correcting  codes  generated  by  the  ideal  I  =<  x±  xg  > 
in  was  creating  a  matrix  that  represented  the  code.  While  a  straightforward  task  for  small 
values,  the  process  was  considerably  speed  up  by  the  following  Sage  code  that  produced  the 
matrix.  The  program  is  described  and  then  examples  are  provided. 


sage:  def  matrix_code_generator(j) : 

II  II  II 


...  A  function  designed  to  create  the  abbreviated  sub-matrix  of  a 

matrix  for  the  code  corresponding  to  the  coordinate  hyperplane  arrangement. 
The  sub-matrix  corresponds  exclusively  to  polynomials  consisting  of  one 
variable  raised  to  a  power.  The  input  is  j,  the  upper  bound  on  the 
polynomials  used  in  generation. 

EXAMPLES:  See  below  for  j=l,...,10 
INPUT :  j 

...  OUTPUT:  factored  determinant  of  the  abbreviated  sub-matrix 


sage 


w=var ( ’  w’ ) 
row=  [0]  *j 

for  k  in  range(j): 

row [k] =w~k 
CodeMatrix= [ [0] * j ] * j 
for  1  in  range(j): 
newrow= [0] * j 
for  k  in  range(j): 

newrow [k] =w~k 
for  m  in  range(j): 

newrow [m] =r ow [m] *  * ( 1+ 1 ) 

CodeMatrix [1] =newrow 
CodeMatrixlist= [0] * ( j  **2) 
for  n  in  range(j): 

for  p  in  range(j): 

CodeMatrixlist [n+ j *p] =CodeMatrix [p] [n] 
MS=MatrixSpace(PolynomialRing(ZZ,w) , j , j) 
Code_matrix=MS .matrix (CodeMatrixlist) 
determinant=det (Code_matrix) 
f actored_det=determinant . factor () 
return  factored_det 
def  matrix_code_echelon(j ) : 

II  II  II 


...  A  function  designed  to  create  the  abbreviated  sub-matrix  of  a 

matrix  for  the  code  corresponding  to  the  coordinate  hyperplane  arrangement. 
The  sub-matrix  corresponds  exclusively  to  polynomials  consisting  of  one 
variable  raised  to  a  power.  The  input  is  j,  the  upper  bound  on  the 
polynomials  used  in  generation. 

...  EXAMPLES:  See  below  for  j=l,...,5 
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INPUT :  j 

...  OUTPUT:  echelon  form  of  the  abbreviated  sub-matrix 

II  II  II 

w=var ( ’  w’ ) 
row=  [0]  *j 

. . .  for  k  in  range (j) : 

row [k] =w~k 
CodeMatrix= [ [0] * j ] * j 
...  for  1  in  range (j): 

newrow= [0] * j 

...  for  k  in  range(j): 

...  newrow  [k] =w~k 

. . .  for  m  in  range(j) : 

...  newrow [m] =row [m] ** (1+1) 

...  CodeMatrix [1] =newrow 

CodeMatrixlist= [0] * ( j  **2) 

. . .  for  n  in  range (j) : 

...  for  p  in  range(j): 

CodeMatrixlist [n+j*p]=CodeMatrix[p] [n] 
MS=MatrixSpace(PolynomialRing(ZZ,w) , j , j) 
Code_matrix=MS .matrix (CodeMatrixlist) 
echelon=Code_matrix . echelon_f orm() 
return  echelon 

sage:  matrix_code_generator (1) 

1 

sage:  matrix_code_echelon(l) 

[1] 

sage:  matrix_code_generator (2) 

(w  -  1)  *  w 

sage:  matrix_code_echelon(2) 

[  1  w] 

[  0  w~2  -  w] 

sage:  matrix_code_generator (3) 

(w  +  1)  *  (w  -  1)~3  *  w~4 
sage:  matrix_code_echelon(3) 

[1  w  w~2] 

[0  w~2  -  w  w~4  -  w~2] 

[0  0  w~6  -  w~5  -  w~4  +  w~3] 

sage:  matrix_code_generator (4) 

(w  +  1)~2  *  (w  -  1)~6  *  w~10  *  (w~2  +  w  +  1) 
sage:  matrix_code_echelon(4) 

[1  w  w~2  w~3] 

[0  w~2-w  w~4-w~2  w~6  -  w~3] 

[0  0  w~6-w~5-w~4+w~3  w~9-w~7-w~6+w~4] 

[00  0  w~12-w~ll-w~10+w~8  +  w~7-w~6] 

sage:  matrix_code_generator (5) 

(w  +  1)~4  *  (w  -  1)~10  *  w~20  *  (w~2  +  1)  *  (w~2  +  w  +  1)~2 
sage:  matrix_code_echelon(5) 
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[1  w  w~2  w~3  w~4] 

[0  w~2-w  w~4-w~2  w~6-w~3  w~8-w~4] 

[0  0  w~6-w~5-w~4+w~3  w~9-w~7-w~6+w~4  w~12-w~9-w"8+w~5] 

[00  0  w~12-w~ll-w~10+w~8+w~7-w~6  w~16-w~14-w~13-w~12+w~ll+w~10+w~9-w~7] 

[000  0  w~20-w~19-w~18+2*w~15-w~12-w~ll+w~10] 

sage:  matrix_code_generator (6) 

(w  +  1)~6  *  (w  -  1)~15  *  w~35  *  (w~2  +  1)~2  *  (w~2  +  w  +  1)~3 

*  (w~4  +  w~3  +  w~2  +  w  +  1) 
sage:  matrix_code_generator (7) 

(w  +  1)~9  *  (w  -  1)~21  *  w~56  *  (w~2  -  w  +  1)  *  (w~2  +  1)~3  *  (w~2  +  w  +  1)~5 

*  (w~4  +  w~3  +  w~2  +  w  +  1)~2 
sage:  matrix_code_generator (8) 

(w  +  1)~12  *  (w  -  1)~28  *  w~84  *  (w~2  -  w  +  1)~2  *  (w~2  +  1)~4 

*  (w~2  +  w  +  1)~7  *  (w~4  +  w~3  +  w~2  +  w  +  1)~3 

*  (w~6  +  w~5  +  w~4  +  w~3  +  w~2  +  w  +  1) 

sage:  matrix_code_generator (9) 

(w  +  1)~16  *  (w  -  1)~36  *  w~120  *  (w~2  -  w  +  1)"3  *  (w~2  +  1)"6 

*  (w~2  +  w  +  1)~9  *  (w~4  +  1)  *  (w~4  +  w~3  +  w~2  +  w  +  1)~4 

*  (w~6  +  w~5  +  w~4  +  w~3  +  w~2  +  w  +  1)~2 

sage:  matrix_code_generator (10) 

(w  +  1)~20  *  (w  -  1)~45  *  w~165  *  (w~2  -  w  +  1)~4  *  (w~2  +  1)"8 

*  (w~2  +  w  +  1)~12  *  (w~4  +  1)~2  *  (w~4  +  w~3  +  w~2  +  w  +  1)“5 

*  (w~6  +  w~3  +  1)  *  (w~6  +  w~5  +  w~4  +  w~3  +  w~2  +  w  +  1)"3 

8.3.  Binary  Skeleton  Codes.  When  working  on  the  j  =  1  case  for  codes  developed  from 
the  simplicial  complex  A(£.  h),  to  prove  the  minimum  distance  formula,  observations  where 
made  concerning  what  the  Hamming  weight  of  the  sum  of  s  row  vectors  became.  Some  Sage 
scripts,  listed  below,  facilitated  the  easy  calculation  of  the  Hamming  weight  of  this  sum. 

def  F(l,h,s) : 
sum  =  0 

for  a  in  ranged  ,h+2)  : 

for  t  in  range (1 , s+1) : 

sum=sum  +  (-2)**(t-l)  *  (binomial (s,t))  *  (binomial (1-t , a-t) ) 

return  sum 


9.  Conclusion 

While  much  has  been  achieved  in  the  work  of  this  project,  room  for  further  research  is  still 
present.  There  are  many  other  classes  of  subspace  arrangement  codes  and  simplicial  com¬ 
plex  codes  where  we  know  neither  the  dimension  nor  the  minimum  distance.  Nonetheless,  by 
varying  the  parameters  £,  h,  and  j  of  the  skeletal  codes  K(£,h,j),  we  are  able  to  produce 
many  useful  codes  with  various  values  of  length,  dimension,  and  minimum  distance.  In  addi¬ 
tion,  it  is  possible  that  the  minimum  distance  and  dimension  formulas  developed  through  this 
project  might  prove  useful  in  understanding  theoretical  problems  in  subspace  arrangements 
and  simplicial  complexes.  Furthermore,  we  examined  a  few  examples  of  our  skeletal  codes  in 
the  McEliece  cryptosystem.  It  has  not  yet  been  determined  (though  it  seems  likely)  that  there 
are  nice  decoding  algorithms  for  our  codes,  thereby  making  them  well-suited  for  the  McEliece 
cryptosystem.  However,  more  work  is  needed  to  provide  a  larger  class  of  codes  with  the  same 
parameters  to  make  the  system  practical.  This  project  was  exciting,  especially  since  it  has 
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provided,  through  open  problems  such  as  Conjecture  6.30,  many  intriguing  future  research 
topics. 
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